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ABSTRACT 

Tiie infrared absorption line shapes involving 
one phonon process have been obtained in cubic crystals 
containing substitutional impurities using thermo dynamic* 
Green's functions and Kubo formalism. With the aid of 
group theoretic arguments it is shown, that only the 
totally symmetric or s-type one dimensional irreducible 
representation of the point group of the crystal contri- 
butes to the impurity vibration and hence, to the optical 
absorption. Conditions for occurence of inband resonant 
modes are investigated. An estimate is made of the 
change of nearest neighbour force constant in the case 

c p 7 

of Li and Li in KBr, based on the experimental data of 
Sievers and Takeno. A formula for the integrated absor- 
ption over all the modes is derived. 

The expressions for the mean square displace- 
ment and velocity of the defect are obtained. The former 

*57 

wa^5 used to compute the Debye-Waller factors for Fe^‘ in 
platinum, palladium and copper and using the available 
experimental data an estimate of the change of force 
constant is made in each case within the framework of 


the models studied 
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The study of the effect of imperfections on the 
physical properties of crystalline solids is an important 
branch of solid state physics. All crystals in nature 
contain some defects or disorders. Presence of defects 
changes the physical properties of crystals in various 
ways. Point defects such as substitutional impurities, 
vacancies and interstitials etc. change the vibrational 
thermodynamic properties of solids by altering the freq- 
uency spectrum. The mechanical properties such as comp- 
ressibility etc. are known to change due to the presence 
of higher dimensional defects such as dislocations (one. 
dimensional defect), stacking faults (two dimensional 
defect) etc. This dissertation contains the results of 
a theoretical investigation of the effect of point defects 
on lattice infrared absorption and Mossbauer effect in 
cubic crystals. 

In crystalline solids in the equilibrium confi- 
guration the atoms are arranged in a regular array called 
the * lattice*. The electrostatic interaction between the 
ions and the electrons is mainly responsible for the 
binding energy of the solid. At finite temperatures the 
ions vibrate about their equilibrium positions* Owing 
to the smallness of the ratio between the electronic and 
nuclear masses, the velocities of the ions will be much 
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smaller than the velocities of the electrons. It is thus 

possible to assume that the electrons follow the nuclear 

motion adiabatically* i.e. the electrons at any time move 

as if the ions were fixed in their instantaneous positions. 

t 

This adiabatic or Bom-Oppenheimer approximation implies 
a potential V(R) which can describe the lattice vibrations. 


lor small vibrations the potential energy 7(R) can 
be expanded in a Taylor's series in terms of the displace- 
ments .^U ^ of the ions from their equilibrium positions, 
lor a Bravais crystal containing N atoms we have 

7(E) = T(E 0 ) + + *2 A a p(J.J) 

u a (i) u p(i) + •• * 


where V ( R 0 ) is the equilibrium potential energy of the 
crystal and 


A a (i) - [ 2Z®. 

j 


i • 

-S = R 0 



r ^ 2 ^cr) . 

Su a (i) ^Up(i).Mo 


(l*2a) 

U#3»> 


r = r 0 + n 




lor harmonic solids the terms up to quadratic in U are reta* 
ined. Higher order terms characterise the anharmonieity* 
In equilibrium configuration the net force acting on any 
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atom is zero, hence the Hamiltonian of the crystal in the 
harmonic approximation is written as 

j p 2 ( i) , 

H = 2j wr- * * /•/ A a g (i,J) <1>3) 

/<* 2 ' j.d, 

j.fi 

where M. is the mass of the ji-th atom. The harmonic force 
constants A a p(J.>.i) reflect the physical symmetries of the 
lattice. The general relationships among Idle force COnSt- 

O^A 

ants are discussed in the standard hooks on the subject? 

In harmonic approximation ‘the system of coupled 
oscillators can he solved exactly with the aid of normal 
co-ordinate transformation which transforms the system of 
coupled oscillators into a set of independent oscillators 
called normal modes, each having its characteristic freq- 
uency of vibration. The translational symmetry of the cry- 
stal leads to a normal mode being in the form of a travel 
lling plane wave, characterized by a wave vector k and 
frequency CO- (k) , j being the branch index. The motion of 
any atom in the lattice is due to the superposition of 
these normal modes. Since the total number of atoms in 
the crystal is very large, the frequencies form quasi-con— 
tinuous bands. This description of solids in the harmonic 

approximation explains many of the properties, such as the 

$ 

low temperature specific heat (Debye T s law), infrared 

, 

reflectivity of ionic crystals, etc*. 
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atom is zero, hence the Hamiltonian of the crystal in the 
harmonic approximation is written as 


j.d. 


2M . 


AttpCifD u a (i> u p(i) 

J,cc 


J ,( 3 

where M. is the mass of the j-th atom, The harmonic force 
JL 

constants A a ^ (i»i) reflect the physical symmetries of the 

lattice. The general relationships among the force const— 

2—4- 

ants are discussed in the standard boohs on the subject* 


In harmonic approximation ' the system of coupled 
oscillators can be solved exactly with the aid of normal 
co-ordinate transformation which . transforms the system of 
coupled oscillators into a set of independent oscillators 
called normal modes, each having its characteristic freq- 
uency of vibration. The translational symmetry of the cry- 
stal leads to a normal mode being in the form of a travel 

* 

lling plane wave, characterized by a wave vector k and 
frequency (k) , j being the branch index. The motion of 
any atom in the lattice is due to the superposition of 
these normal modes. Since the total number of atoms in 
the crystal is very large, the frequencies form quasi-oon- 
tinuous bands. This description of solids in the harmonic 
approximation explains many of the properties, such as the 

ft 

low temperature specific heat (Debye T* law), infrared 

i ° 

reflectivity of ionic crystals, etc*. 
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In the language of quantum theory the lattice vibra- 
tions can he treated as a field and quantization of such, 
a field re suit s in what are known as 'phonons** The lattice 
vibrational quanta or phpnons have no widths in the harmonic 
approximation and hence have infinitely long life time, hut 
inclusion of anharmonic effects brings about interaction 
among the phonons leading to a finite life time which plays 
a vital role in lattice thermal conductivity, thermal expa- 
nsion, shape of fundamental absorption line in ionic crystals 
and a host of other properties. She phonon description of 
lattice vibrations is closely related to the classical normal 
mode analysis, because the independent normal modes get 
coupled in presence of the non-linear terms in the crystal 
potential energy. One switches over to the phonon picture 
by expressing the normal modes in terms of annihilation and 
creation operators. However, for higher order processes the 
phonon description is more suitable as there exist convenient 

tools such as many-particle Green's functions and diagraram- 

6 

atic techniques of the perturbation theory for the many- 
body problem. Besides the phonon-phonon interaction arising 
out of anharmonic ity, interaction of phonons with other ele- 
mentary excitations also adds to their life time. 

She interaction among the phonons is also brought 
about by the defects which are present in most of the solids. 
The pioneering work in this field was that of Lifshitz^*"' 1 ’^ 
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who had evaluated the perturbed frequencies and normal modes 

in an imperfect crystal and discussed the other physical 

11 12 ■ 

properties of such a crystal* Montroll and Potts' 9 gave 
a method of attach for this problem based on ihe use of 
Green's function of difference operators, and they also app- 
lied it to some idealized models of lattices with impurities* 
Montroll 1 ^ et al , Tanaka and Yamahuzi^ Mahan ty, Maradudin 
and Weiss'^' 1 *’ and the Russian school lead by Bifshitz^ have 
studied the effect of impurities on vibrational thermodynamic 

properties which are additive functions of the normal mode 

18 

frequencies. Rubin has studied the detailed dynamics of 
an impurity atom in a lattice. 

Besides the change in the normal mode frequencies, the 
localized nature of the impurity perturbation depending on 
its strength (mass of the impurity and its interaction with 
the host) gives rise to two interesting features. Sew modes 
of vibration, are often generated having frequencies above 
the maximum frequency of the unperturbed crystal, accompanied 
by a certain change in the density of states within the bands. 
These modes are called 'localized modes* as the displacement 
amplitude of such a mode falls off with increasing distance 
from the impurity site. In harmonic solids such a high 
frequency mode cannot decay into the band and has a 6— function 
line shape. However, anharmonicity^ and impurity-impurity 
interaction * leads to a finite life time for this mode 
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and hence a weil defined shape to the line* Secondly, under 
suitable conditions an impurity does not produce a local mode, 
but a significant redistribution of levels within the bands, 
occurs in such a way that in the neighbourhood of a certain 
frequency a marked enhancement in the density of states takes 
place, such a peaking up in the density of states is known as 
a 'resonance’ or ’quasi-localized’ state* Due to the presence 
of a large number of states in the vicinity this mode decays 
with a finite life time even in the harmonic approximation* 
Localized modes can also occur in the forbidden regions 
between the acoustic and optic branches of the frequency 
spectra in ionic crystals* These modes are called the gap 
modes* 


The existence of localized modes in solids was known 

from the work of Lifshitz, though the possibility of occurence 

of inband resonant modes was over-looked in the earlier stages 

of the work. Brout and Visscherf^ Kagan and I o sile vskii^ and 
24 

Takeno independently had pointed out this possibility theo- 
retically. Brout and Yisscher were able to calculate the 
frequency of this low lying mode caused by a heavy isotopic 
impurity (impurity host interaction remains unchanged) in 
Debye model as 



where ZlM = M - M = mass of the impurity - mass of the host atom 
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^ = Debye cut off frequency* 

The width of this node was calculated in detail* 

Koster and Slater^*^* had considered the effect of 
an impurity atom which produces a different potential at the 
impurity site, on the electron energy spectrum of a metal* 

The problem of scattering of spin wares due to impurities 
in ferromagnets was discussed by Wolfram and Callaway ^ ’ ^ 

OQ 

and Takeno independently. There exists a great deal of 
similarity in the problems of impurity affecting the electr- 
ons, the phonons and the magnons in solids which makes it 
possible to treat them on the same mathematical footing* 

This aspect has been discussed by Izyumovf^ 

The problem of lattice dynamics of imperfect crystals 
inrolres the study of the abore aspects both theoretically 
and experimentally. Both the localised and resonant modes 
were first predicted theoretically* With the arailability 
of experimental techniques in the measurement of low tempe- 
rature thermal conductivity, specific heat, elastic 
inelastic scattering of neutrons, lattice infrared absorption 
and Mossbauer effect one can go into the details of defect 
vibrations. Here we briefly outline the various experimental 
' methods for studying the impurity activated modes. 

i 

In case of insulators the thermal resistance is pro- 
duced due to the scattering-^"*^ of lattice vibrations or 
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phonons. Although the impurities are knoijn to scatter phonons, 
the main contribution to thermal resistance at high tempera- 
tures comes from the phonon-phonon scattering due to the an- 
harmonic terms in the crystal potential energy. However, at 
very low temperatures the impurity contribution sometimes 
becomes much more pronounced, so as to be observed distinctly. 
Pohl^ first observed a dip in the low temperature region of 
the thermal conductivity curve when plotted against tempera- 
ture for potassium chloride crystal containing a small concen- 
tration of potassium nitrite, The usual Raleigh type of 
scattering' of phonons at low frequencies could not expl- 
ain the result and Pohl reproduced the dip by assuming a 
resonant form for the inverse relaxation time. Subsequently 
the effect was confirmed by a series of measurements of the 
thermal conductivity of KOI doped with KI, UaCl, and CaClg by 
Walker and Pohli^ 

The occurence of Brout-Visscher resonance diie to a 
heavy mass defect or softening of impurity-host interaction 
led to the conclusion that the resonant scattering of phonons 
at low frequency which in temperature scale corresponds to 
such a low temperature was responsible for the thermal condu- 
ctivity dip. This was the first experimental evidence for 
the occurence of inband resonant mode in solids. In case of 
a molecular defect like and CN“ the resonance origin is 

due to the internal degrees of freedom of the molecular 
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vibrations whose frequency lies embedded in the phonon 
spectrum. The molecular impurity problem was discussed 
theoretically by TJagner-" who introduced molecular 
Green* s functions characterising the new degrees of freedom 
and solved the resulting set of coupled equations for 
scattering and localized modes. 

Since the defect changes drastically the frequency 
spectrum* the properties such as specific heat which directly 
depends on it should show up the impurity model® In fact, 
the low frequency resonant mode which produced a marked 

enhancement in the specific heat at low temperature was 

41 42 

detected by Lehman et al and Panova and Samoilov in Mg 

containing Pb and Cd as impurities. 

In perfect crystals, the inelastic coherrent and 
incoherrent scattering of neutrons involving one phonon 
process directly reveals the dispersion (frequency wave 
vector dependence) relation and frequency spectrum respect- 
ively in cubic crystals. The experimental detection of the 
low frequency resonant mode was made by M/ller and Mackintosh 

in chromium mixed with tungsten and by Svensson, Brockhouse 
44 

and Rowe in copper doped with gold. The first group of- 
experimenteers also observed a temperature independent shift 
( non-anharmonio ) of the resonant mode. Such a shift which 

9 

arises due to impurity- impurity interaction was calculated 

45 

by Elliott and Maradudin ^ by considering many isotopic 
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impurities distributed, at random in the crystal. The result 
could not agree quantitatively with the above experiments, 
although there was a qualitative agreement. Recently Behera 
and Deo^ extended the theory by including the change in 
impurity host interaction and using diagrammatic method in a 
rather phenomenological way. To our knowledge the discre- 
pancy between the theoretical calculation and result obtained 

from the experiment still exists. Elliott and Taylor^ and 
48 

Taylor discussed the solution of Dyson equation for many 

4.Q 

impurities in co-ordinate space. Recently Takeno obtained 
a formal solution in case of an extended defects (including 
change of impurity host force constants). Although these 
are improvements over the previous theories, the results are 
far from satisfactory for higher concentrations of impurities. 

The existence of localized modes were observed as peaks 
in the energy distribution of neutrons scattered inelastically 
and incoherrently in palladium containing nickel as light 
impurities by Mozea* Qtnes and Myers, ^ 

Interaction of electromagnetic radiation with matter 
is by far one of the oldest and fundamental branches in 
physics. In crystalline solids this interaction yields 
information about its structure and throws light on mechanisms 
of various physical processes. The existence of a non-vani- 
shing dipole moment in harmonic crystals corresponding to 
the transverse optic mode of zero wave-vector-' (cQ^q) is 



1m own to cause Reststrahln phenomena in ionic crystals. As 
regards lattice vibrations whose frequencies extend from zero 
to a certain maximum value which lies in. the far infrared 
region of electromagnetic spectrum, the optical absorption 
experiments provide perhaps the most direct method available 
at present, although the experimental difficulties of far 
infrared spectroscopy are substantial, 

Schaefer^ 1 observed for the first time in I960 a 
high frequency lattice mode in the infrared in alkali halide 
crystals containing hydride ions as impurities, This is known 
as U-center and is caused by the appearance of an infrared 
active local mode outside the maximum frequency of the crystal. 

9 

Schaefer also observed sidebands on both sides of the main 
U-center band. Hie frequency of the U-center was calculated 
bj Wallis and Haradudin, and by Takeno - ^ et al , using a 
nearest neighbour force model for the host crystal and obtai- 
ned a good agreement with the experiment in the mass defect 
approximation. However, when long range electrostatic forces 
are taken into account the agreement becomes poor. More 

C A 

realistic calculation by Jaswal and Montgomery using regid 
ion and deformation dipole models, and by Fieschi, Hardelli, 
and Terzi using a shell model showed that a 50 /f reduct- 
ion of nearest neighbour force constant is required to fit the 
data. Mirlin and Reshina^ and Fritz'*®*'^ have studied the 
details of the wing structure and the temperature dependence 
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of the m„-in hand. Pritz observed a strong temperature dep- 
endent peak vith half width decreasing as T 2 at room tempe- 
ratures and he ascribed this and the appearance of the side 
-ands Lo tho interaction of the local mode with the band 

modes through the anharmonicity of the crystal potential 

fin — fii 

energy. Timusk and Klein have developed a theory to 

explain the side band structure in terms of the anharmonie 

~ 62 
coupling of H ion to its nearest neighbours. Nguyen has 

performed a similar calculation including both the anharmonie 

and second order electric dipole moment interaction. The 

later result showed that the dominant factor arises due to 

anharmonicity . Experimental and theoretical investigations 

6*5 

were also carried out by Mitra and Brada ^ whose results were 
similar to that of earlier observations* 

fiA fi £> fifi 

Balkanski and Nazarewicz, ** 3 and Ingress et al 
observed infrared active local modes and band modes, acti- 
vated b; : , the impurities like boron and phosphorus in silicon. 

67 

Montgomery and Hardy 1 studied the isotopic variations of 
dispersion frequency in mixed alkali halides. 

As regards the inband resonant modes in ionic cry- 
stals are concerned, the first observations were those of 
68 

Sievers in potassium halides containing silver ions as 
impurities. These resonances lie in the low frequency region 
of the acoustic band ( 100cm to 10cm ) and arise due to 
the heavy mass of the silver ion. He in collaboration with 
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Takeno could detect a very low lying sharp impurity mode 

for very light impurities like Li * ' in KBr. In all these 

cases the absorption frequency decreases slightly and half- 

70 

width increases with the rise of temperature. .Their obs- 
ervation gave an impetus in understanding the configuration 
of small ions like Li in alkali halide crystals and this 

will be discussed in Chapter IV. Localized modes inside Idle 

71 

forbidden region or the gap modes were also observed by 
Sievers and his co-workers in KE due to Cl - and Br”. 

The existence of molecular impurity centers in 

alkali halides were known from the thermal conductivity 
^"5 

measurements. However, thermal conductivity results alone 

cannot furnish all the informations regarding the internal 

degrees of freedom like free rotations, hindered rotations 

and tunneling etc.. The molecular modes are infrared active 

and show up in the absorption experiments. The energy levels 

of molecules and transition from hindered rotation to free 

rotation with increase of temperature were studied by Seward, 

72—74 

Harayanamurti and Pohl. A somewhat clearer picture of 

molecular vibrations in solids emerges from the combined 
study of both these properties. 

Infrared measurements were also carried out in rare 

73 

gas crystals by Jones and Woodfine who measured the one 
phonon absorption spectra of solid argon containing nearly 
2 of krypton or xenon (neutral impunities) • Very recently 
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Klein arid Macdonald^ 6 detected, the isotopic absorption spec- 
trum in ilc.01 and LiB crystals. The natural Gl^ f ^ isotopes 
acting as simple mass defects result in absorption in the 
acouscic band and yield information about the phonon frequen- 
cies at critical points. This is the first observation of 
isotope induced band mode absorption in alkali halides to our 
knowledge . 


The theoretical work in this field is mainly confi- 
ned to the isotopic impurity case or the so called mass defect 
approximation (MDA) • In this approximation an expression for 

absorption co-efficient due to charged Impurities in otherwise 

77 

uncharged host was obtained by Dawber and Elliott. In ionic 

crystals the notables are the work of Hardy, Gunther^ and 

iq 

Llaradudin. Liaradudin has given a general formula for the 

absorption co-efficient using the propagator technique. 

Recently the case of an extended impurity (involving a change 

of mass as well as nearest neighbour force constant) was 

80 81 

discussed by Patnaik and Mahan ty, Benedek and Nardelli 

Q p 

and Takeno. 

Before discussing any further regarding the defect- 
activated absorption we will briefly discuss the mechanism of 
interaction between the crystal and the radiation field. In 
the process of absorption, emission or scattering of electro- 
magnetic radiation, the quantum mechanical description involves 
two states of the crystal namely the initial and the final 
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‘state. Here we assume that both these states do not take 
into account the motion of electron inside the ions or atoms. 
In otherwords we exclude from our consideration the elect- 
ronic transitions. Hence the infrared absorption is concerned 
with the process where the initial state is the crystal in 
the ground state and an incident photon, whereas the final 
state is an excited state of the crystal with one or more 
.phonons, This process can take place if the crystal possesses 
dipole moment. Like the potential energy, the dipole moment 
depends on the nuclear co-ordinates and hence can be expanded 
in terms of their displacements from the equilibrium position 
in a Taylor's series as 



Z e ji yi>Vi> + 

i>P 


2> Y 


(1.5) 


where M„(R) = a- cartesian component of the dipole moment of 
the crystal, 


and 



J>M„(R) 

[ — ] ■ 

• DupQ) R = R 0 


9 



(r) 

— a ] 

(i)'£u (J) -R = R c 


(1.6a) 


(1.6b) 


^ e^ j and J satisfy conditions 8 similar to these of 

force constants, 3?or cubic crystals M a (R) is isotropic and 
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K° which, represents a permanent dipole moment does not Occur, 

(X 

hence we can write 


M a (R) = Zj e j u a (i) +H e ^ /U a ^ u a^ + ••• 


22 


( 1 . 7 ) 


The co-efficient e. represents the ’effective charge’ of 
"til 

the j, atom and the tern containing it is the linear moment. 
Ihe co-efficients e.. which are involved in the second order 

_L.T_ 

moment, determine the effective charge of j atom which is 

*tll 

produced by the displacement of jj_' atom. In ionic crystals 
tne first order dipolemoment exists whereas the effective 
charge is zero for homopolar crystals such as those having 
diamond structure. The second and higher order moments exist 
for both polar and non polar crystals. 


lor polar crystals the existence of the linear moment 
leads to one phonon absorption process only at one frequency 
T - vi ' Ul a S- function shape. This is the well known fund- 
amental absorption line, sometimes called the Reststrahl 
absorption. Although there is no one phonon absorption in 
crystals having diamond structure, they do exhibit a weak 
well-defined absorption spectrum involving two .or more phonons 
due to higher order moments. It had been first proposed by 
Lax and Bur stein that -the presence of charged defects in 
homopolar crystals will make the whole range of lattice modes 
optically active. In ionic crystals the defect destroys the 
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translational symmetry and relaxes the selection rules. 
However, it is not necessary to introduce charged substitu- 
tional impurities, as any aton becomes polarised by distor- 
tions induced by the vibrations of the neighbouring atoms 
due to their unequal electron affinity. The induced first 
order dipole moment interacts with the electromagnetic 
radiation to give optical absorption, the experimental evi- 
dence of such absorption spectra were observed by Jones and 

7 15 

Wood-fine in rare gas crystals as mentioned earlier and 

Q/ 

by Newman and Willis ^ in silicon containing carbon as 
impurities. 

Theoretical calculation of effective charges were 

ojr 0 

made by Szigeti in homopolar crystals and by Hartman 

87 

and Elliott ' in rare gas crystals. Before closing the 

formal discussion on the photon-phonon interaction, it may be 

pointed ov/c that induced enelastic scattering of light by 

phonons (Hainan effect) takes place in ionic crystals due to 

the presence of a defect which changes the electronic polar- 
88 

izability of the crystal in first order. 

89 

After the discovery of Mossbauer effect i.e. the 
recoil free emission or absorption of y-rays in crystals, it 
became a very useful tool for the study of defect dynamics. 
There are three dynamical quantities of interest in connect- 
ion with the resonant absorption or emission of y-rays 
in crystals. The 'zero phonon' process in which the crystal 
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remains in the same state after the absorption (or emission) 
of photon is known as Mossbauer effect. The quantum mechan- 
ical probability for occurence of such process is given by 
the Debye-Waller factor which depends on the mean square 
value of the displacement of the emitting nucleus. The sec- 
ond quantity of interest is the shift (temperature dependent) 
of the position of the zero phonon line due to an increase 
in the vibrational energy of the crystal, because of decrease 
in mass of the emitting nucleus as a consequence of the 
relativistic equivalence of mass and energy* This is known 
as second order 'Doppler shift' and is proportional to mean 
square velocity. The third quantity relates to 'one phonon 
process’ in which the y-xa.jB photon is absorbed with a change 
of energy of the crystal by one quantum of vibration. This 
gives information about the frequency distribution function 
of the ho-': crystal when the resonant nucleus does not differ 
from the rcom :.t replaces (trie isotopic impurity case). The 
resonant mode or local mode if present will show up in one 
phonon process. The literature in this field is very vast* 

The works of visscher,^ 0 Llaradudin,^ 1 and Lehman and DeWames^ 
deal with the defect dynamics associated with the Mossbauer 
effect. Most of the earlier papers on the experimental and 
theoretical status are discussed by Draunfelder."^ 

The brief review given here is not intended to be 
exhaustive. Detailed reviews of the experimental and theore- 
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tical situation are available 
articles . 19 - 9 *' 99 


in a number of review 


The present thesis is mainly oonoerned with the 
theoretical investigations of impurity activated infrared 
absorption due to a charged substitutional impurity in an 
otherwise uncharged host, and by an impurity in ionic cry- 
stal taking into account the change of mass as well as the 
nearest neighbour force constant in some solvable models of 
cubic crystals. Expressions for the mean square velocity and 
displacement of the impurity are obtained, and the results 
compared with the available experimental data, A brief 
summary of the thesis is as followsj 

In Chapter II we will discuss the mathematical 
formalism that is to be used. The double time temperature 
dependent Green's functions are discussed in general without 
assuming any specific form of the impurity potential. 
Expressions for the displacement and velocity co— relations 
functions are obtained in terms of the above Green's funct- 
ions, le use the Kubo formula to study the response of the 
system to an external time dependent perturbation like that 
due to incident electromagnetic radiation and a formula for 
optical absorption (Power absorption) is derived# A sum rule 
for the total integrated absorption including the band and 
localized modes is obtained and applied to the case of charged 
impurities and impurities in ionic crystals. 

* t 
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Chapter III deals with the absorption line shapes 
of charged inpurities in simple cubic, body centered cubic 
and face-centered cubic crystals. She symmetry properties 
of the defect lattice is taken into accou_it to decompose 
the resonant denominator , into various factors representing 
the different irreducible representations of the point group** 
By group theoretical argument it is proved that only the 
totally symmetric or s-wave mode is infrared active. lor 
understanding the low frequency mode a Deybe model is emplo- 
yed to obtain expression for the frequency of resonant mode 
and its half width. 

In Chapter IV infrared absorption in ionic crystals 
is considered. The model of ionic lattice is taken as a 
simple cubic lattice having alternatively positive and nega- 
tive charges associated with the two different masses. Using 
the M transformation, which makes it a simple cubic lattice 
with an effective mass, the absorption line shape formula 
is discussed. The portioned forms of the perturbation and 
Green’s function matrix are employed to sum over all the 
lattice sites. The absorption due to impurity is being 
separated from that of the host crystals* The results are 
compared with the experimental data of Sievers and Takeno. 

In estimation for the change of nearest neighbour interaction 
is made. 


The mean square displacement and velocity are derived 
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for an inpurity in Chapter V. An estination of inpurity host 
nearest neighbour force constant is male by fitting the theory 
to the experimental measurements of Debye Waller factor for 
le in Pt, Pd and Cu lattices. In all these models a good 
fit to the experiment is obtained for a wide range of tempe- 
ratures. A series expansion formula is derived for the mean 
square displacement which give also the correct fit. 

In the concluding chapter the results of the 
investigation are summarised. The limitations of the theory 
is also discussed. 
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CHAPTER II 

GREEN’S FUNCTION FORMALISM AND EXPRESSION FOR POWER ABSORPTION 
2.1 General Formalism 

It has been already pointed out that the entire solid 
from the point cl view of lattice vibrations in harmonic 
approximation can be considered as an assembly of phonons. 
Hence the study of the vibrational properties of solids is 
essentially a many-body problem. The defects bring about 
interaction among the phonons and scattering takes place. 

Many formalisms were developed to treat such problems. The 
usual perturbation theory is not applicable to the study of 
localized modes. The defect problem is exactly solvable in 
harmonic approximation by making use of the classical equa- 
tions of motion of the lattice and the Green’s function 
technique devised by Lif shitz^® , Montroll and Potts • The 
same problem can also be solved by quantizing the lattice 
vibrational field and then using the Green’s function tech- 
nique as prescribed by quantum field theory. However, both 
the methods are equivalent for the assembly of harmonic 
oscillators. 

Here we use the method of the double time thermo- 
dynamic Green’s function which has already been applied 
successfully and extensively for a variety of interacting 
systems. But the Hamiltonian for the lattice vibrations 
and hence the Green’s functions are expressed in terms of 
the space co-ordinates rather than in the momentum space* 

In this form, the displacement and velocity correlation of 



any two atoms can be expressed directly in terms of their 
Green's functions. 

In literature there exist many reviews and textg on 
these Green's function formalism. We shall follow olosely 
the review written hy Zubarev*" ^ There are three different 
types of Green's functions defined in terms of the average 
value of the I product of the operators, namely the causal, 
the advanced, and the retarded. We use only the retarded 
Green's function defined as 

; = <A(t), B(t )> = - B(t')]> 

, ( 2 . 1 ) 

where A(t) and B(t ) are any two Heisenberg operators i*e# 

A(t) = e iHt A(o)e" iHt (2,2) 

(a system of units in which h = 1) , H is the Hamiltonian of 
the system. The square bracket denotes the commutator 

[A,B] = AB - BA (2.5) 

and 6?(t-t ) is the Heaviside step function, 

t 

(t-t"") * 0 for t > t' 

t 

= 1 for t < t' (2.,4) 

The angular bracket denotes the ensumble average, i,e. for 
any operator A, 

<^A^> = Tr e'^A/Tr e~^ H (2*5) 

where p =(kgT) , kg being the Bolt man constant and T is 

the temperature of the system. The equation of motion obeyed 
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by the operator A is 

i = [A,H] ( 2*6) 

Differentiating equation( 2.1) with, respect to t we obtain the 
equation of motion for the Greeks function as, 

i H = &(t-t ) <[A(t) r B(t )]> + ^[A(t),H], B(t - (2.7) 

The second term in the above equation is another 
Green's function which is generally of higher order (involving! 
more number of operators) than the original one. The usual 
way is to write the equation of motion of that Green's funct- 
ion which in its turns leads to a hierarchy of coupled equat- 
ions. To solve this set of never ending equations one has to 
decouple the Green's function at a certain stage* There are 
no set rules which will decide this temination, but one has 
to make suitable approximation taking into consideration the 
physical situation. But in case of the defects lattice 
problem, the second Green's function is of the same order as 
that of the original one and its equation of motion yields 
the first Green's function. Hence one has to solve the 
resulting coupled equations simultaneously. 


Let G(^’) be the Fourier transform of the Green's 


function G(t-t ), then 


G(t-t ) = 

- 00 


and 


(t-t ) 0(to) auJ( 




G(ccj) = ( e 1 ^ ) d(t-t' ) . 


(2*8a) 


(2*8b) 


-05 
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She physical quantities of interest in many problems 
are the time correlation functions defined as 

F BA (t,t' ) = <B(t ) A(t)^> . (2.9a) 

*&(*,*') = ^ A(t) ' (2.9b) 

Equal time correlation functions are called the auto-corre- 
lation functions and they give the average value of the prod- 
uct of the operators. One can express the correlation funct- 
ions in terms of the spectral density function J (<^0 ) 


+ CO 


4 f f 

F BA (t-t ) = dujJ(<o) e“ ic0 

(2.10a) 

-* -t-cc ~ - , 

^AB ( t~t ) = Ucoe^' J(co) e~ ±u) K 

(2.10b) _ 

J 


G(cC +ie^ - G(w-ie) = -i (e^-l) J(<t>), 

(2.11) 


and 0 is implied. 

Ihe abo”o expression relating the Green's function 
and the spectral .. r.sity function is important. Once the 
retarded Green's function is known, the spectral density 
function can he obtained and then taking the Fourier transform 
one can get the expressions for various correlation functions 
at all times and at all temperatures. 

§ 2.2 Linear Response Sheory And Kabo Formula 

jf 

„ , . 101.102 

Kubo has given an elegant formulation of the problem 

of reaction of an external perturbation* on a quantum mechanical 
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system while discussing transport phenomena like the elect- 
rical and thermal conductivity, Peltier effect etc* in terms 
of the time correlation functions of the respective currents* 
The theory is cased on the assumption that the response of 
the system is linear to the external perturbation (i*e« elec- 
tric field in case of electrical conductivity etc*). The 
higher or the nonlinear terms are neglected in the above 
analysis, but they are important in the theory of phenomena 
which arise due to non-linearity* 

let us consider a system with time independent Hamil- 
tonian H. A perturbation h| which depends on time explicitly 
has been switched on at t = - a» , when the system was in ther- 
mal equilibrium* How the system is described by the Hamilto- 
nian 

E= H + * (2.12) 

The ensemble average of expectation value of any operator A 
is given by 

< A > = Tr( f(t)A) - (2.13) 

where / P(t) is the density matrix representing the ensemble 
and satisfy the following equation 

= [H, jt(t)] = [H+H*, f(t)] . (2.14) 

The above equation is to be solved under the initial condition 
(-00) = f> = e - pH /Ir e^ H , (2.15) 


to the first order of perturbation, i.e* 
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jP(t) = /-" L ^p(’k) 
Zi^Kt) = 0 for t = - £o 


(2.16a) 

(2.16b) 


Restricting the solution upto the linear term in the 
external perturbation and by making use of the initial condi* 
tion the change in the density matrix is obtained as 

Afl t) = - i j > [H* ,/>] . ( 2 ' 1 7) 

Row the average value of any operator A is given by 


A = Tr (( p+£fK t))A) 


= A 0 +i f^L( t), h| (t )J> dt . 


(2.18) 


This expression gives the response of the system to the exter*- 
nal perturbation in terms of the retarded Green’s function* 

^ 2.5 Electrical Conductivity 

The crystal interacts with the electric field of the 
electromagnetic radiation through its dipole moment, let the 
electric field E(t) of radiation be uniform throughout the 


entire space and vary sinusoidally with time, then 
E(t) = E cos .Ot 


(2.19) 


The corresponding interaction Hamiltonian is 


Hj(t) = -M.E(t) = "2^ e jj_ E*u(j,t) cos^t e £ 


( 2 . 20 ) 


where e is a small positive infinitesimal quantity indicating 
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the fact that the interaction was switched on at t = -to. 


Under the interaction there arises an electric current 
and if one considers a unit volume of the crystal, the expect- 
ation value of the current density operator "becomes 

. /X. 


where 


with 



ft 

- =j \fj a (t), H*(t' )^> dt'. 

-O0 1 

(2.21) 

Hj(-fc) 

= Ej(t) cost^t e £ ^ , 

(2.22a) 

J a (t) 

=2 e i ^ a,t) ' 

J 

(2.22h) 

Hj(t) 

= - 2 e j E a u a ( ^'^ • 

(2.22c) 


Integrating "by parts and using equation (2.22a) this reduces 
to 


<J a OO) = - Ee| <[J a (t), H J (t)]> 


i (i«J+e)t 


co -3-e 


> f 00 


- *)> 


.(i 1 ^ +~)t 
iw+ e 


dt 


( 2 . 23 ) 


cj 


The first term involves a commutator at equal time and is 
evaluated with the help of the usual relation. 


[ u a Ci5* )3 = stt 5 ap 6 14 ' j 

J 


(2* 24a) 


and 


Since 


[u a (j?}> )] - E^ a (3)» ^U)3 = 0 . (2.241s) 


S 


I (t) = -£ e i \ ^ (2 - 25) 


J* 

equation (2.23) is rewritten as 
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<J 0 (t)> = Re £ c ajS to) Bp e (iw+ e) 


(2.26) 


where the element of thr conductivity tensor is given hy 

+ iC 


a ap (io) = ' £ H? 


cep 



^(o), J p (t 


- 00 


(iyJ+ e) t , 

® dt 

iw + e 


(2.27) 


The first tern in the above expression corresponds to 
the electrical conductivity of a system of free charges. The 
-conductivity tensor can also be expressed in terms of the time 
correlation of currents. However, it is expressed here in 
terms of the retarded Green’s function directly, lor conve- 
nience we define the velocity Green’s function which is rela- 
ted to current Green’s function by a multiplying factor and is 
similar to that of Elliott and Tayloa? - , 0 ^ 

<^u a U,0), . (2-28) 


In terms of its Eourier transform the conductivity tensor is 
expressed as 





ejSw 

A 


0 1 


^p(i.iW) . 

co + xe 

(2.29) 


^ 2.4 Green’s Functions And Correlation lunotions Of The lattice 

In order to calculate the lattice electrical conductivity 
arid electric current density etc. the knowledge of the lattice 
Green’s functions is required. The expression for the Green’s 
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£ unction is obtained "by writing down the equation of motion 
of the Green’s function as defined hy the equation(2,28) . i.e. 


1 h - ’ 






(2.30) 


Differentiating once again 


jaL i> 

dt 2 


= - 6 (-t) 




2 ff 0PyUfS* t ) 

n,/ 


\ 6 (£>i) 

( Vj ) 4 


(2.31) 


Here we have used the commutation relations as given by equat- 
ions^. 24) and the lattice Hamiltonian as defined in Ghapter I. 
The above equation can be solved by matrix inversion method. 

In matrix form, this and the defining equation(2.28) are writ- 
ten respectively as 


and 


’* , . A A A 

•=- P = - 6 (-t) D - PD , 

dt 2 



( 2 . 32 ) 


(2,33). 


The cap denotes the matrix. U(t) is the 3N dimensional vector • 
whose (a,m) element denotes the a-eomponent of the velocity of 
the m^k atom. 


where 


D = 


Jo-4 £o jJo-4 


and hence ^o 


0 


(2.34a) 

(2.34fc) 
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The matrices i, ft 0 . A and *° are at* dimensions. A 
is the harmonio force constant matrix and K is tbs diagonsi 
oass matrix. The superscript denotes tbe quantities of «• 
unperturbed crystal. The eigenvalues of D° are the squares 

of the normal mode frequencies ^<£) of the unperturbed 

• fA a A Vitr q wave Vector h SUld "branch 
lattice and they are specified hy a a - 

index j. The eigenvalue equation is written as 

Yj V £3 5 = ^ 3 ^ 9a ^ ' ■ * 2 * 355 

ffi/9 

, f m / l whi ch are normalized plane waves 

The eigenvectors wuen axe ^ 

u «v 0n +innq satisfv the following 

travelling along x, y, a directions, satisiy 

orthonormality and olosure relatione. 


£ %!■$ ,p d% ) = 6 S' 8 33 - ’ 
£», <* 

E W = S 2£ ‘«P - 
*,/ 


(2.36a) 


(2.36b) 


(2.37) 


The explioit form for <P a (£§) ls 

. (.8) = — e^S^e (M) - 

Vfca' • fE 

where -fe(hj) V are the polarization vectors. Their orthenor- 

mality anl olosure properties follow from equations(2.36). 1« 
Brails crystals e(kj) are real. The normal mode frequencies 

satisfy the condition, 

Wj (-k) =cu,(£> • C2 - 38) 
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The plane wave character of the lattice vibrations 

is due to the translational invariance of the crystal.* In 

the presence of defects this is destroyed and the vibrations 

depart from their plane wave nature. There are two distinct 

ways of looking at the problem for its solution. The usual 

method is to express the vibrations, the Green* s functions 

etc. in terms of the unperturbed one's whose properties are 

well known. The alternative way is to look for a suitable 

transformation which can diagonalize the perturbed dynamical 

matrix. The eigenvalues of this matrix are the squares of the , 

normal mode frequencies of the perturbed crystal and hence^ 

cannot be labelled by the propagation vector k and branch 

index 3 . In stead, they can be labelled by an index s whioh 

runs from 1 to 3H. In practioe such a transformation is ■ 

difficult to obtain explicitly. However, one can obtain a > 

formal solution of the problem which can be used later to j 

' ! 

derive certain other useful relationships. The eigenvalue i 

I 

equation of the perturbed lattice is \ 

\ 

23 Cfi.5) 9, < I ) = co B < ( 1 ) (2.39) ! 

fi.fi f 

- ' 

where -£c<p s j are the new normal mode frequencies. The eigen- ! 
vectors satisfy the following orthonormality and closure j 

relations. 

Z, < ( 1 > < ( &> " 6 B3 ' C 2 .«a) 

rn, <x. . 

S < ( • ) ( S > - & aa s ap • (*•«»> 

s , 
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It is convenient to use the first method for most of 
Our calculations. The Fourier transform of equation (2.32) 
yields 

* 

F (to) = — + F ( -O ) I) > (2.41) 

271 

This is written in a more convenient manner as 

F ('O ) = - — [ i - u) z G ] (2.42) 

2tz 

A 

where I = is the 317 dimensional unit matrix and 

G = (L!oj 2 - Af 1 . (2.43) 

✓v 

G is the Green’s function matrix whose properties are dis- 
cussed by Montroll and Potts. This Green’s function was 
widely used to study the dynamical properties of pure crystals 

as well as crystals containing substitutional defects. Recen— 

104 

tly this technique has been applied to molecular lattices 

and other impurity clusters which occur naturally as disorders 

105 

involving a line or plane of defects. The surface defects 

which gives an additional mode known as ’surface mode’ has 

106 

also been discussed in terms of this Green’s function. The 
unperturbed Green’s function is given by 

G° = (Mo/ - A) 0 " 1 (2.44) 

Any particular element of this Green’s function is expressed 
as 

( S , a; * ± *> - <vy-fc 

3 
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JL 

3N(M M )£ 

Hi - j . 


s— , e g (kj) e s (^j) 

4-' “* - fcj ± le 


i^*(S(m)-R(n)) ^ 

(2.45) 


The summation over k is over the first Brillouin zone* 


Now one can obtain easily the velocity and displacement 
Green* s functions in terms of G • Taking the Fourier trans- 
form of the equation( 2*28) and equation( 2*42) we have 



Hence -the displacement Green* s function 



(2*46a) 


(2*46b) 


The time correlation functions of the Green's functions are 
obtained with the help of equations (2 *10) and equation(2.11) 


u a (m y o)^> = liffl 


J 

e-^ o -oo 


+ (O 

f da ) 

Ce^-l) 


[ & a p(S »£» +ie) - G aR (m, n^J 3 - is)] (2.47a) 




. A3 

„i*’t , . 

e dw> 


<5? <s» *) u a <S'°4> = 

E -*0 - 00 

[G a p(m, nj^ 2 +ie) - (^(mjnjcj 8 - ie)] 


ap 


(2.47b) 


It is clear from these expressions that in the harmonic 
approximation the Green's functions do not depend on tempera- 
ture and hence any properties which are directly related to 
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Green’s function 11600120 temperature independent. But the 
correlation functions do depend on temperatures. Ihe general 
formula for correlation functions are valid for all times and 

k i • 

for all temperatures. By making m = n F a == p and t = o, the 
respective auto correlation functions are obtained and they 
express the mean square velocity and displacement of that 
particular atom. Ihese two quantities are important in conn- 
ection with Mossbauer effect where the resonant atom is gen- 
erally a substitutional impurity in the crystal. 


$ 


2.5 Power Absorption 


In presence of electromagnetic radiation the energy is 
continually transferred to the crystal due to interaction 
through the electric dipole moment. She strength of the abs- 
orption depends on the frequency of the exciting field and 
the nature of the phonon spectrum of the crystal. One of the 
convenient Wo,ys of studying the response of the system to the 
external radiation is to measure the power lost by the field. 
Ihe expression for power absorption B(to) per unit volume ad 
a function of exciting frequenoy is obtained by taking the 
time average of the scalar product of the eleotric field with 
the current density i.e. 


B<W)'-[Kt). <£(*)> Average 


003 ^ Es 2>e <6 ° ) E P Average 

oC p 

( 2 . 48 ) 



..36 


■With tiie help of equations(2.29) and (2,42) it oan be expre- 
ssed in terms of the Green's function as 

B(«0) = ^ ^ e m e n E cc E p 111 G a p(S*£:»c<> S “ ie ) • (2.49) 

n 

80 

Shis is an exact expression for the absorption line- 
shape function for a harmonic crystal when the absorption is 
due to one phonon prooess only. No specific model of the 
lattice is assumed in deriving this formula. She applicabil- 
ity of this formula depends on the evaluation of the Green's 
function and they are model dependent. In almost all the cases 
(except Debye model) the Green's functions are to be evaluated 
numerically with the help of high speed computors. 

As has been pointed out earlier the harmonic Green's 
functions are independent of temperature, and hence the infra- 
red absorption comes out to be temperature independent. 
However, in reality i.e. measurement of absorption co-effici- 
ent temperature comes into play a major role. A temperature 
factor can be incorporated in the theory if one considers 
higher. order phonon processes arising out of the anharmonic 
terms in the potential energy or the non-linear terms in the 
electric moment or both together. But difficulties arise in 
solving the problem of defect together with enharmonic ity or 
higher order dipole moment, as it yields a number of coupled 
equations involving higher order Green's functions. 


The above expression which involves only one phonon 
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process can be obtained by solving classically the Newton’s 
equations for the displacements of the lattice* In the 
presence of radiation, the charged impurities in otherwise 
uncharged host or ions in polar crystals undergo forced 
vibrations with the same frequency as that of the radiation, 
but there is a phase shiftf’°^between the driving force and the 
lattice vibrations. The expression for power absorption as a 
function of the exciting frequency is obtained by talcing the 
time average of the scalar product of velocity of the charged 
atom with the driving force and this depends on the phase 
angle. Since the phase shift is related to the Green’s func- 
tion, the ultimate result comes out to be same as that of 
equation(2.49) . 

6 2.6 Sum Rule 

71 

Dawber and Elliott have derived an expression for the 
total integrated absorption over both the band and localized 
modes. The contributions of both types of modes were calcu- 
lated separately. For localized modes, they have assumed a 
normalized shape function which takes into account the broad- 
ening of the level due to anharmonic effects. However, one 
can derive such a result for a general defect which involves 
a change of force constant besides the change of mass in the 

I 

harmonic approximation without assuming any particular shape 
for ihe localized mode. 

a 

T 

The total absorption B is defined as 
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B* = 


f 


CO 


B(co) duJ 



e m e n E a I m G a p(m,nj^ s -ie)dcu . 

( 2 ,' 50 ) 


As has Been pointed out earlier that the perturbed dynamical 
matrix can he diagonalized by a unitary transformation, so 
that we can express the perturbed Green’s function directly as 


:_i r 

W7 ""V 16 


(2-51) 


Henoe 


B 1 = 


e e 

7C m n 


1 L, 

f Jjr$ 




V P £ 


(|)<y sco/ 1 -^) 
( 2 . 52 ) 


Now the 6-function integral is carried out easily. She 
closure relation given by(equation( 2.40b) helps us to write 

B *“ f • <2.53) 

fr f o{ ~ 


This formula is quite general irrespective of any 
specific model of the crystal lattioe, henoe can be applied 
to the case of charged impurity as well as impurity in ionic 
crystals. No parameter characterising the interaction bet- 
ween the atoms occur in this expression. Only parameters 
which determines the total absorption, depends on the 
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impurity mass and its charge, 

(i) Charged Impurities: 

For a dilute concentration of charged impurities 
(U-j. per unit volume) whose masses are in otherwise unchar- 
ged host crystal, the total absorption is given by 

e 8 E s II t 


B 


I 


% 


= T 


M' 


(2.54) 


(ii) Impurities in Ionic Crystals: 


We crnsider the ionic crystal as consisting of 
two types of masses (M^ and Mg) with charges -e and +e at 
odd and even sites respectively. If impurities with 
charge +e and mass m[ are substituted at the odd sites the 
total absorption comes out to be 

(2.55) 

66 

It has been pointed out by Angress et al. that 

the effective charge of an impurity can be estimated from the 

total area of the absorption curve provided the concentration 

08 

of the impurities is known. Maradudin has also obtained a 
similar expression for the total absorption. 
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CHAPTER III 

ABSORPTION LINE SHAPE DDE TO CHARGED IMPURITIES 

^ 3*1 Lattice Model And Perturbation 

For a perfect harmonic lattice the amplitude of 
vibrations are plane waves. It is well known that in the 
presence of defect due to the loss of translational symmetry, 
the vibrations depart from their plane wave characters. 

Since the perturbation is local one, the modes are strongly 
perturbed near the vicinity of the defect, and hence these 
modes can be studied best by using as a probe some property 
of the defect itself. In case of crystals which do not 
possess a first order electriG moment the introduction of an 
impurity with non zero effective charge leads to absorption 
of photon by creation of one phonon over the entire range 
of lattice modes. 

A charged substitutional impurity in a crystal ’ 
excited by monochromatic radiation is one of the simplest 
example of a forced oscillator coupled to a large system. 

An oscillator coupled to a system with many degrees of free- 
dom and continuous spectrum behaves like a damped oscillator, 
the damping being due to the continual transfer of energy 
from the oscillator to other modes of the system. The absor- 
ption line shape of a charged substitutional impurity would 
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•therefore, be intimately related to the structure of the 

frequency spectrum of the crystal with the impurity. In 

reality the .optical absorption occurs through one and 

multi phonon processes, but the one phonon process is of 

primary importance, the multi phonon process leading to a 

background against which the one phonon absorption is 

thrown into relief. It is known that absorption spectra 

108 

of type I diamond arises due to nitrogen impurities. 

Recently there have been considerable interests 

in the field of phonon scattering by defects, specifically 

due to the presence of quasi-localised mode in the acoustic 

oo 1 OQ 

band which mani-fests in resonant scattering. Callaway, * y 

lakeno,^ Wagner,^ 8 Krumhansl,^‘ I '*' >, ' J " 1 ' 1 Thoma and Ludwig"*"^ 

79 

and Gunther J have considered the details of the problem 
with particular interest in the resonant scattering. 

Krumhansl and Mathew considered the problem of one dimensional 

10^5 

lattice in a rather complete way. Elliott and laylor 

used the double time thermodynamic Green’s function to solve 

90 

the isotope defect. Visscher has discussed a different 

aspect of the defect vibration which is of importance in the 

Mossbauer effect. Yussouff and Mahanty^ 1 ^ * and Kleil'^*^ 3 '^ 

have given general theories of scattering of phonons in cubic 

109 

crystals-. Group theoretic method was applied by Callaway, 
to separate out the scattering amplitude into its various 
symmetry parts characterised by different irreducible 
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representations of tlie point group of the lattice and it 
was followed up by the other worlrers. ill these studies 
lead in a natural way to the exaot evaluation of the Green’s 
functions of the perturbed lattice when the impurity atom 
involves both a ehange of mass as well as the nearest neigh- 
bours force constant. 

We consider here the case of one defect in the 

crystal and its effect on the vibrational properties. In 

reality the crystal contains a number of impurities, but 

in many physical situations where the impurity concentrar- 

tion is very low one obtains a good approximation for many 

impurities in the crystal by multiplying the effect produced 

by a single impurity by the total number of impurities* 

This is the case for example in impurity induced infrared 

absorption in ionic and homopolar crystals where the impur- * 

18 5 

ity concentration rarely exceed 10 per cm . In this 
process one neglects the impurity-impurity interaction which 
gives rise to effects in higher order of concentration. 

Knowledge of perturbed Green’s function is indis- 
pensable in calculating various properties of impure crystals. 
This depends on the structure of the crystal as well as the 
interacting forces among its constituent atoms. One cannot 
proceed to calculate the detailed behaviour of a crystal 
without considering any specific model. Any realistic model 
becomes almost intractable due to heavy numerical computations 
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which sometimes obscure the physical situation to a consi- 
derable ext ent • 

We consider monatomic cubic lattices (simple cubic, 
body-centered cubic and face-centered cubic) and investigate 
the cases where the polarization does not enter into the 
picture. The lattice consists of N atoms each of them 
having mass M. The harmonic interactions are limited upto 
the nearest neighbours in such a way that the three cartesian 
components of the displacement vector become independent of 
each other. In this so called scalar model the phonons are 
not polarized. Such unrealistic lattice models are chosen 
for mathematical simplicity in analysis, however many of the 
features of optical absorption by impurity in this model are 
of sufficient generality to be applicable to more realistic 
models. For the study of low frequency impurity modes whioh 
are being discussed here, such a simple but mathematically 
tractable model is particularly useful because the low 
frequency behaviour of phonons is not so sensitive to the 
crystal model. 

The lattice sites are denoted by a vector R given 

by 

Rg - la e-j_ + ma eg + na e^ , 

Rg = (1+ Tp)a e^ + ?ja ©g + (n+ *?j)a * 

Rp = £(l+n)a e^+jKl+nOa e 2 +Hm+n)a e^ . 


(3.1a) 

(3.1b) 

(3*lo) 
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Here a is the lattice constant and (l,m,n) are any integers 
(zero included). (o-j_,S c ,e^) are unit vectors along x,y and 
z directions respectively. She subscript S,B and ]? denote 
respectively the simple cubic, body-centered cubic and face- 
centered cubic lattices. The time independent equations of 
motion of the lattice in the matrix form is 


(M^ 2 - A) U = 0. (3.2) 

where the force constant matrix A has the following structure 
A s (l,m,n/l,m',n) = y[6 6-^- 6^ 6^ - 6^ 6^- 

^ 6 l+lil + 6 1-1,1^ 

“ 6 11‘ 6 nn' (6 m+l,m' + 6 m-l,m /) 

“ 6 11 6 mm ^ 6 n+l,n + 6 n-l,n^» (3.3a) 

A B (l,m,n/l,m,n,) = y[3 6-^ 5^ 6^- 6^ 6^ 

^m,m+l + ^ia,m-l^ 

“ 6 nn^ 6 l,l+l 6 m,mr-l + 6 i',l-l 6 m>+l 5 ' 

^ll^m,m-l ^h,n+ 1 + ^m^m+l ^n,n-l^ 

“ (6 1,1+1 6 m,m-l 6 n,n+l + 6 1 / ,1-1 6 m',m+l 


(3.3b) 
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and 


^(l,m,n/l,Li,n) = y' 126^/ 6^ 

“ S 11 £ nn ( 6 m,m+l +6 m',in-l)" 6 mm 6 nn ^ 6 l',l+l + 6 1,1-1^ 
- ^n y } n-l + ^m^m-l £ n,n+l^ 

^mm^3^L+l £ n',n-l + £ 1,1-1 £ n,n+l^ 


6 r,w( 6 n / ui _ 1 + 5-,' -i -i 6 / _ .,)]. 

iui 1,1+1 m,m-l 1,1—1 2 n,m+l 


(3.3c) 


The solutions of equation (3.2) are the plane waves given by 
U(R) = p e 1 ^'^ . (3.4) 


Substituting equation(3.4) in equation(3.2) and making use 
of the equations (3. 3) the frequency of wave-vector relations 
are obtained as 


<o| (& = 
ioi <*> - 


(£> = 


^ [ 3 - cos k x - cos k 2 - cos k^] ,(3.5a) 
^ [ 1 - cosj^k^. ooss-kg.cosik^] , (3.5b) 

|p- [ 3 - cosi-k-L.cosikg-cosi-kg.cosl'k^ 

- cosl-k^.cosi^] . (3.5c) 


It is convenient to assume that the substitutional 
impurity whose mass is M+AM is substituted at the origin. 
The nearest neighbours are coupled to it by force constant 
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Y+^Y* So the perturbation is limited to the impurity site 
and its nearest neighbours. The model of a simple cubic 
lattice with impurity is shown in Pig. 1. In presence of 
this perturbation the equation of motion becomes 

(Mu; 2 - A) U = [ (M«.' 2 - A) 0 + P] U = 0. (5.6) 

A 

The elements of the P matrix are given by 


P(0,0) = M<o 2 - nny , 

. V = - ay . 

P(0,R £ ) “ p (V°) =AV • (5.7) 

where n is the number of nearest neighbours. It has value 
6 for simple cubic, 8 for body-centered cubic, and 12 for 
faoe-centered cubic respectively. 


The equation relating the perturbed and unperturbed 
Green’s function follows from e quat ions ( 2.43) and (3.6) as 


G- 



- G°PG 


(3.8) 


lor many impurities the perturbation matrix P has several 
blocks and the equation(3.8) is the Pyson’s equation, For 
many impurities the exact solution is difficult and one 
solves the problem by iteration. The formal solution of the 
above equation is 

5 - (i + n)- 1 ^ i 


(3.9) 



SIMPLE CUBIC LATTICE MODEL WITh IMPURITY 



. FI 
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For ore impurity the perturbation matrix has only a small 
number of non-vanishing elements constituting a sub matrix 
the dimensionality of which is determined by the extent of 
the perturbation. This helps us to write P matrix in the 
following partitioned form 


P = 



0 

0 

0 



( 3 . 10 ) 


where p is the non-vanishing block in the perturbation matrix. 
Its dimensionality is 7,9,13 for simple cubic, body-centered 
cubic, and face-centered cubic lattices respectively. The 
unperturbed C-reen’s function matrix can be partitioned in a , 
similar manner i.e. 


V 


s 

i 

i 

1 e '21 

s 12 

s 13 \ 


§22 

§23 * * ♦ \ 

(3.11) 

\ g 31 

S 32 

A { 

g 33 ... j 


\ * 

• 

m • *> • j 

/ 



A Q 

where each of the sub-matrices of G has the same dimen- 
sionality as that of p. With the help of these two equations 
(I + G°P) takes the form as 
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r A '*■ 

I+GP 


■' 1 + SP 


Son P 


3 21 
’31 


P 


0 

A 

I 

/i 

0 


c 


( 3 . 12 ) 


. / 


The inverse of this matrix is easy to work out, because 

A ^ ^ 

it involves the inverse matrix of (I + gp) whose dimension- 
ality is small. Henoe the perturbed Green's function 
becomes 


f (f+gp r 1 

! 

A 

0 

- \ 

0 

\ 

/ A A 
/« s 12 

f 

A 

g 13 

! -g 2 iP(i + gp) _1 

A 

I 

\ 

0 ... 

/ ®21 g 22 

g 23 

6 = ! -gjxPCl+gp)" 1 

A 

0 

i ...! 

®31 g 32 

g 33 




I 


(3.13) 


This is the general form of the perturbed Green's 
function matrix. Any element of it which is required for 
the specific calculation can be obtained by considering the 
appropriate block of this matrix. In fact, as far as the 
optical absorption due to charged impurities is concerned 
one needs some definite element which makes it simpler to 
select the block. 

^ 3.2 Absorption Line Shape 


Since the impurity having charge 'e* at the origin 
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we have 

= c- 5 . n . (5.14) 

4- 4i.» — 

Then the absorption line shape formula given by equation(2.49) 
reduces to 

B(t-c?) =i s 2 E 2 6 o5h 0 (C,0, ^. 2 -is) (3.15) 

It is clear now that only the (0,0) element of the perturbed 
Green’s function enters into the calculation. This can be 
obtained from the 1st block of the matrix given by equation 
(3.13). 

0(0,0, v.' 2 -ie) = [(I+gp)” 1 ^]^ . (3.16) 

This element 3an be calculated by using the explicit 
form of p and g matrices and then inverting the matrix 
(I+gp) directly. However, this inversion can be achieved in 
a more elegant way by taking into account the symmetry of 
the system. This method not only saves the labour, but 
throws some more light on the nature of the vibrations of 
the impurity which would have been otherwise obscure in the 
direct calculations. 

Although the defect destroys the translational 
symmetry, the point group symmetry of the crystal about the 
impurity site as the reference point is still retained. This 
enables one to use the displacements of all the affected sites 
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in forming a reducible representation P of the point group 
of the lattice, which can be decomposed in terms of irredu- 

-y i > 

cible representations j 'of the point group as 

f - £ a, r' W • (3.17) 

2’ 

where a , are zero or integers. The matrices p and g have 
the symmetry of the point group of the lattice and hence, 

A 

one can find out a unitary matrix V which bloc£ diagonalizes 

p, g, and any function of them, namely (I+gp). Such a 

matrix can be found by group theoretic procedure. The number 

117 

of such blocks, the dimensionality of each block etc., are 

' i 

j}/ 

known from the equation(3.17) . If J has the dimensionality 
n^ , then in block diagonal form there will be n^ identical 
a p x a^ blocks. The cubic crystals whose point group symm- 
etry is Oj, the reducible representation p' is decomposed into 

its various irreducible representations in the usual notat- 

. 118 
ion as 


Is ■ = 

A lg 


F lu + E g * 

(3.18a) 

Ib - 2 

A lg 

+ 

E lu + P 2g + A 2u * 

(3.18b) 

= 2 

A lg 

■f 

F lu + 3 g + F 2g + E 2u * 

(3.18c) 


This shows the various symmetry vibrations possible in the 

A 

cubic lattices. The V matrix is constructed with the help 
of projection operator and is given in the appendix I. The 

* >\ S' 

rows and columns of the p, g, and V matrices are labelled 



..52 


in same manner. For example in case of simple cubic lattice 
the rows end cnlurms are labelled as (000), (100), (010), 
(001), (-10C), (0-10) (Go-1). 

The block diagonal form of (I+gp)~ for simple 
cubic lattice is given by 



0 T" 1 / (3.19) 

•*_1 

where the various blocks are as follows. T^ is a two 
dimensional matrix belonging to s_wave (JL ) part. Its 

-*-o 

explicit form is 

J-l = n / 1 ^ S 1 -' 6 ^(Soa-gj.) 

1 D S \ -T6!2c'*g 1 (^ - ^C) 1+ ^ +(1- Wg os 

(3.20a) 




and 

Tg 1 = = T" 1 = ~ , (I lu ) p- wave parts. (3.20b) 

Dg 

T" 1 = Tg 1 = ^ , (Eg) d- wave parts. (3.20c) 

D g 

where Det(I+gp) = l I + gp = Dg(l|) 8 (Dg) 2 . 


(3.20d) 
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The nomenclature s, p, d (used here as superscripts) 
corresponding to the irreducible representations A^ g , 3? lu , 
and E are used to describe the partial wave amplitudes of 

O 

vibration. They have acquired such names in analogy with the 
symmetry of s, p, d wave functions in a cubic crystal. The 
explicit forms of the factors of the determinant are 


Dg = 1 + SX + *1 (l + dg) _ »J[L fljjd* + 

f^T ( 1+ T[)]g os * (5.21a) 
= 1 - M (g os - gjj.) . (3.21b) 

d| =1 - Ay (g QS + g 2 - 2g l;L ). (3.21c) 


The four different elements of the Green's function 
used above are given below. 


g os = G°(000, 000; CO 2 - ie) 

gl = G°(000, + 100 i<o 9 -ie) 

(Any permutations of the indices) 

S 11 = G °(± 10 °* 0 ± 10»'to 2 -ie) 

S 2 = G °<± 100 » + 100;OJ 2 -ie) 


(3.22a) 

(3.22b) 


(3.22c) 

(3.22d) 


Evaluation of these Green's functions are possible by nume- 
rical methods which are well known. 
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In the above calculation contains all the four 
different elements of the Green’s functions, "but the existence 
of some identity relations among these elements simplifies 
this to its present form. Ihese relations can he deduced 
from the definition of the Green’s function matrix for the 
perfect lattice i.e# 

(M CO 2 - A) G 0 = I • (5.23) 

Or taking element on both sides, we have 


£ 


- h 1>a 




(3.24) 


Now by making use of the equation(3.3a) and equations (3. 22) 
the following relations are obtained 


V-luf - 6y) g QS + 6y g x = 1 , (3.25a) 

(lliLC z - 6y) g 1 + y (g os + 4g 11 + g 2 ) = 0 . 

(3.25b) 

From equations (3. 20) it is clear that only the 
s-wave part of the determinant involves both the parameters 
( ^|| , ) of the perturbation. Ihe p and d parts depend 

only on the change of force constant. In case of isotopic 
impurity only this mode of vibration is possible. With 
this decomposition of the perturbed Green’s function, its 
(0,0) element is obtained as 
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& S (0»0»4) 2 -ie) 


V* A A ~l A 

[V Q T J -V ( 


t g] 


S &J 0,0 



^(1- 


MU) 2 


s os 


1 £i 
QY Y 




(3.26) 


The interesting point to he mentioned here is that, 
the (0,0) element contains only the s-wave block correspond- 
ing to the totally symmetric one dimensional irreducible 
representation of the point group of the lattice, hence this 
mode becomes infrared active. This statement is valid not 
only for cubic lattices, but for all lattices -which are 

monatomic and have the symmetry of any one of the point 
107 

groups. This is proved in appendix I. The pysical rea- 
soning lies in the fact that except the s-eave mode the 
impurity does not take part in other modes of vibration 
although its nearest neighbours do move. This is because the 
p-like (three fold degenerate) and d-like (two fold degene- 
rate) modes have nodes at the impurity site. . Since the 
impurity atom is charged the dipole moment contribution 
comes only from the s-wave mode. 


Since the columns of the V matrix are normalized 
vectors which are obtained from a general vector character- 
izing the displacements of the affected sites by projection 
operator, they represent the specific symmetry mode vibrat- 
ions of the affected sites. The various symmetry modes 
belonging to the different irreducible representat ions .-given 
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■by equation(3.18a) are shown schematically in Fig- 2. 


Since the other two modes are not infrared active 
they cannot he studied by optical absorption experiments. In 
this respect the impurity and its nearest neighbour vibrati- 
ons resemble a molecular vibration problem. The (0,0) ele- 
ments of the perturbed Green’s functions for bGC and fee 
lattices can be obtained in a similar manner, ill these 
above discussions and arguments hold in those cases and we 
write them as 


u/-ie) 


(l+(l- §^ 2 ) ^ jg 0E + |y 


-JC,(%27a) 


& p(0,0i c0 2 -ie) 



01 

Y 


(3* 27b) 


where 



= 1+ 01 + £L 
Y Y 


Mj . dM\ 

' x+ M ; 8y 




- 4jj (1 + 5X)] gaB , (3.28a) 

D| - 1 + « + (l + ~ E ^U +<a 5> TW 

(1 + «C)]g OT . (3.28b) 


The present forms of equations (3. 28) are obtained by using 
the following relationships among the elements of the Green’s 
function and they are similar to equations (3. 25) 



5 



FIG. 2 
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For too, (Me o z - 8y)g QB + 8y = 1, (3.29a) 

(UiV z - 3y)g 1B + y (g QB + 3g 2B '+ 3g 3B + g 4B ) = 0 • 

(3.29t) 


where 

g 0B = G b (0,0,0, 0; to 2 -ie) (3.30a) 

g iB = 0 B (£»£»£, 0;w 2 -ie) , (3.30h) 

g 2B = 0 B (l,0,0,0,u; 2 -ie) , (3.30o) 

g 3B = G° (l,l,0,0,u; 2 -ie) , (3.30d) 

g 4B = G° (1,1, 1,0, oj 2 -ie) . (3.30e) 

For fco, 

(Mio 2 -12y)g 0F + 12y g ip = 1 , (3.31a) 


(Mcj 2 -12y)g 1B + y(g 0B +4g- LB +2g 2B +4g 3 j,+g 4B )-0. (3.31t) 


where 

g 0 , = Gj (O.O.O.O.^-ie) , 
g^_ B = (.itht 0,_0 i (_fj —ie) , 
g 2B = Op (1»0,0,0; q) Z ~ is) , 

S 3F = & F * 

g 4B = (l,l,0,0,<o 2 -ie) . 


(3.32a) 

(3.32h) 

(3,32o) 

(3.32d) 


(3.32e) 
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The (C,_C) elements in all these three cases are similar 
excop c the structure factor and their Green’s functions, 
hence one can write a single expression for the line shape 
function for power absorption as 


3(~t;) = i e 2 3 2 (l+ 



sQ Im gp 

[Roh s ] 2 + [lmD s ] 2 


(5.33) 


where Re D s and Imh s denote the real and imaginary parts 
respectively of the resonant denominator D . Taking the 
computed values of g Q for various frequencies the power 
absorption can he plotted for the three different crystals. 


3.3 Discussion 

The poles of the perturbed Green’s function deter*- 
nines the now features in the spectrum of the lattice. 
Considering the equation (3*9) we obtain the equation 

Det (I + G°P) = jl + g pi =0- (3.34) 

l 1 

which gives all the characteristic frequencies of the impu- 
rity modes. This determinant is a product of several 
factors as discussed earlier. Ill of them are wholly real 
outside the band and the roots (if any) lying in that region 
axe the localized mode frequencies. Since there is no 
imaginary parts, the life time of those modes are infinite. 
Those local modes belong to different symmetry as they are 
obtained from different partial wave parts of the determinant. 
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Inside the hand these equations are complex, the 
real parts of which may vanish depending on the strength of 
the perturbation. If the imaginary part is small at a fre- 
quency for which the real part vanishes, a resonant mode is 
created. One can cal enlace the life time of such a mode in 
the follov.ing way. 

The density of states for the squarred frequencies 
is given by 


/ r '(^) =7ir 2 suo 8 -«;*) = -jgj 23-r 


2 2 
Sw -co- 3-e 

5 


(5.35) 

Using ^he orthogonality relation (2.40) this reduces t6 

Mu 8 ) = ^ Hr (kH #) , 


= — [M TrG°+rMG(0,0) -M Tr G^I+G 0 ]?)" 1 ^ 0 ] . (3.36) 


- 1 ', 


With some manipulations the expression becomes 





[M TrG°+ — log) .I+G 0 P| ] 
du.> 2 


(3.37) 


where 

MIrS 0 = f- Im g 0 (co z -iE) . (3.3ft)- 

*= density of states for squared frequencies 
for unperturbed lattice. 

The change in squared density of states is 
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&p(co z ) 


Im d 

¥rc * lo § 


i + gpj 


IS Z 


Ntc i 


(3.38b) 


where i takes all the factors of the determinant. These 
additional terms superimpose on the unperturbed density of 
states. The piculiarities arises due to the denominators 
D 1 . Suppose there is a resonance in the i-th channel, then 
the change in the density of states in this mode is 

AP(l .») =1 Im 4 -i- . (3.39) 

u pi Sic [ReD 1 ] 2 + [imD 1 ] 2 

If id the resonant mode frequency i.e. ReD^C^) =0, then 
we can expand ReD 1 in the neighbourhood of as 


/. 

ReD 1 = ReD 1 ^) (io 2 - to\ ) + ... 


(3.40) 


The higher order terms are neglected. So the approximate 
fnrm of equation (3.39) in the neighbourhood oft^ reduces 
to a Lorentzian shape. 


/\p. (af) = V — 1 — - . (3.41a) 

with a width given by 

f / . 2 \ = _ ImD 1 (o; 2 1 # (3.41b) 

li(! ° * [ReD 1 ^ 2 )]^ 

Knowing the resonance frequency one can calculate the width 
of the resonant level. 
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In view of the above discussion the absorption line 
shape function can be expressed in terms of the frequency 
distribution function as 


b(ul> ) = gf ffi 

4M [ReD s ] 2 + [lmD s ] 2 


(3.42) 


where 

f-'U O) = 2a) / C(£o 2 ) 

= frequency distribution function of the 
unperturbed lattice. 

The denominator D s gives the additional structure 
to the phonon density of states. In case of an usual atom 
( ^ = 0, QL =0.) being ionized the absorption spectrum 
directly reveals the frequency distribution function. This 
enables one to determine the one phonon density of states for 
the host crystal if the impurity is really isotopic and 

~ 0, 3ut the line shape function undergoes drastic modi- 
fications when the impurity is accompanied by a change of 
mass or force constant. 

The absorption line shape due to an isotopic impu- 
rity in simple cubic lattice is illustrated in Fig. 3. The 
unperturbed density of states is shown by the curve Fig. 3a. 
For heavy mass defect there appears a pronounced peak in the 
low frequency region which shifts further towards lower edge 
of the band with increase in mass. With the shift the 
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sharpness also inoreases. For very light Isotopes the 
absorption spectrum readjusts itself (Fig. 3o and 3d) with 
localized modes appearing above the maximum frequency of the 
host lattice. For an isotope defect whose mass is slightly 
less than the normal atom a high frequency reso nan t mode 
(Fig. 3b) is produced. In this case the position and shape 
of the absorption line is very susceptible to the perturbat- 
ion parameters. In Fig. 4 the variation of the shape and 
location of the above line is shown with respect to a change 
of force constant. By increasing the impurity host force 
cons bant by 20 '/. the line transforms into a very well defi- 
ned one at the band edge of the high frequency side. Further 
increase in force constant or decrease in mass, shifts it out 
of tho band. A 20 ’/ reduction in the impurity host force 
constant shifts the line towards the low frequency side with 
a broadening cf its width. But further decrease in the force 
constant the shape becomes finer and moves towards the lower 
edge. 

These resonances are due to the peaking ups in the 
density of states as has been shown in equation (3.38b) and 
pointed out first by Brout and Yisscher and Takeno. This 
specific mode also shows up in the thermal conductivity dip. 
Inband resonant modes for heavy mass defects or tremendous 
softening in the force constant are common for all the latt- 
ice s we have investigated. The difference in their positions 
etc. arises due to the structure of the lattices. 
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In case of bcc lattice model 7 /e have considered, 
there exists a logarithmic singularity at the mid point of 
the band (co 2 A) 2 = .5) . The unperturbed density of states 
goes to infinity at that point resulting zero absorption. 
Various absorption curves are shown in Pig. 5 and 6. Pig. 5a 
illustrate the density of states of the unperturbed lattice. 
The usual sharp absorption line arises in the low frequency 
region (Pig. 5b) for heavy mass defect. Por light isotopes 
(Pig. 5c) the spectrum readjusts itself. The effect of 
impurity host interaction is plotted in Pig. 6. The curve 
(Pig. 6a) illustrate the case of an isotope impurity whose 
iH.aos is QC '/. of that of the host atom. There are two peaks 
on both sides of the mid point, but the peak on higher 
frequency side is the stronger one. By readucing the force 
constant by 20 '/ the low frequency side peak (Pig. 6b) grows 
stronger with a simultaneous shift towards the lower edge, 
whereas the high frequency side peak gets diminished* By 
enhancing the force constant by 30 '/ the reverse effect is 
obtained. The high frequency side peak forms a well defined 
line at the band edge whereas the low frequency side absorp- 
tion becomes almost flat (Pig. 6c). 

The absorption line shapes for foe lattices are not 
plotted as they show similar low and high frequency absorption 
lines. The high frequency absorption can also be achieved 
by a pair of sufficiently light impurities as has been shown 
by Tokeno 11 ^ in simple cubic lattice. 
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Since most of the experimental observations of 
impurity absorption modes were detected in the very iLow 
frequency region of the frequency spectrum one can employ 
the isotropic Debye model to evaluate the Green’s function. 
Dor a simple cubic lattice 


He g 


os 


- — [l+i-s^ log "tLt 1 | -I 

^ ^ri 3 


Im g os = & 


3% 1 
2 


M^> 2 


CO 

wj 


U H +1 

COScD, 


m 


= o , 


CO^cO. 


'•/here 


r*i 


<o: 


12v 
I * 


(3.43a) 

(3.43b) 


(3.43c) 


For very low frequency region one can get the analytical 
expression for width. The resonance frequency is obtained 
by putting ReD s = 0, and using equation (3.43a). To lowest 
order in.-j 2 this fields a solution for the resonant mode as 
2 


■Jj, 

CO' 


(1 + ) 

, hM AM , £X _ 2 0I 
5 TS + M Y Y 


(3.44a) 


With the help of this value of the resonant mode frequency 
the width of the mode is given by 


!m > - f • *3 >w - <’•*«» 

0 i 

Takeno and Sievers 6 ^ have calculated these quantities to fit 
the experimental data in case of mixed alkali halides. 


n 
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CHPAT3R IV 
IOUIC CRYSTALS 

.<*■ -4«1 Lattice Model And Green* s Function 

In ttds chapter we consider the response of alkali 
halide crystals (JTaOl type lattice) to electromagnetic 
radiation in the infrared part of the spectrum. It is a 
characteristic property of ionic crystals that they have 
a strong absorption in the infrared associated with the 
motion of charges of opposite sign toward each other. In 
the framework of harmonic theory the existence of a first 
order dipole moment at zero wave vector corresponding to 
the transverse optic mode (called the dispersion oscilla- 
tor) is solely responsible for such a phenomena. Besides 
this strong absorption maximum the experimental observation 
reveals a weak continuous absorption throughout the spectrum 
which arises mainly due to the anharmonicity of the crystal 
potential energy. However in case of mixed alkali halides 
(either anion impurities or cation impurities, sometimes 
both together) the Fourier transform of dipole moment exists 
for all the frequencies making the entire spectrum optically 
active. 

The space lattice of ionic crystal (UaCl structure) 
is face oentered cubic with one cation and one anion with 
each lattice point. On e constructs such a crystal structure 



"by arranging alternately tile two different ions at the 
lattice points of a simple cubic lattice. The features of 
the model studied here are as follows. The lattice consis- 
ting of M ions having masses and Mg with charges e and 
-e are situated alternately at the even and odd sites 
respectively. The interaction is limited to the nearest 
neighbours through equal central and noncentral forces with 
force cons bant equal to y. The charge distribution is 
assumed bo be spherically symmetric and hence the total 
charge of the ion is assumed to be concentrated at the 
center. This assumption implies that the charges move 
rigidly without any distortion. 

The impurity whose mass is and nearest 

neighbour force constant y+/\y is situated at the origin. 

120 

This model is similar to that of Mitani and Takeno. The 
time independent equations of notion of the lattice are 
given by 

(M .u' 2 -6v)x(jj > )+y[x( j^+l, dj.?d 2 

+x(3 1 ,n 2 -i,d 3 )+x(j 1 , j 2 , j 3 +i)+x(j 1 , 3 2 ,3 3 -i)] 

+ H £(i*i) x (i) = °* (4.1a) 

1 

where 

xQ) = displacement of j-th ion, 

1 = (3i>32*3 3 ) ^ 3 ~ 


(4.1b) 
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element of the usual perturbation matrix 
described in Chapter III. Equation( 4 . 1 a) gives two diff- 
erent equations with and M2 for even and odd sites 
respectively* Both types can be combined into a single 
equation by M transformation and the resulting equation 
in the matrix form is 


[(mL> 2 -A) + P ] U = 0. (4.2a) 

with 

MiO 2 = [(M-jW 2 - 6y) (M 2 oj 2 - 6y)]*, (4.2b) 

uQ) = (M -^60 2 - 6y)^xQ), for 3= even, ( 4 . 2 o) 

= (MgU) 2 - 6y) 2 xQ), for 3 = odd. (4. 2d) 


A is the usual force constant matrix for a simple cubic 
lattice. The new perturbation matrix P is modified as 

P(j.,j.') = P PQ,i), if d = e ven » j'= even, (4.3a) 

= “ if d = odd, d /= ° dd (4.3b) 

= if odd » 5 /= even 3c) 

d = even, Y - odd 


where 


p = [ (M 2 ^ s - 6y)/(Mf) z - 6y)] 2 


( 4 . 3 d) 
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Now the diatomic simple cubic lattice can. be 
treated as a monatomic cub^c lattice with effective kinetic 
energy. The other discussion follows the same line as that 
oi simple cubic lattice. The solution of equation(4.2a) 
with P ~ 0 are the plane waves. So the Green’s function 
becomes 


- ±e) = [M(o 2 - ie) - A]“ : 


= 1 JP . (4.4) 

1 "X Moh 2 - Mco 2 (k)-ie 

^ * — 

The asterisk refers to the lattice formed by the M trans- 
formation. The frequency wave vector dispersion relation 
is obtained by solving the equation 

MtJ 2 -Med 2 (k)=[ (M.jtd 2 -6Y) (M2^<- ,2 -6Y)] ^ +2Y(cosk 1 +cosk2+ 

cosk2) = 0 .(4.5) 

The solutions are 

]* 

(4.6) 



There exists two regions for which the frequency 
wave vector relationship is satisfied leaving a gap inbet- 
ween through which no wave can propagate. The upper and 
lower sign Inside the bracket in the above equation give 
the optic and acoustic mode frequencies respectively. Here 
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= 6 t (MJ 1 + ic; 1 ) +< 0 ^ (4.7a) 

"■ Maximum frequency of the lattice i.e. 
top of the optic branch. 

For , 

O 2 0 = (Bottom of the optic branch), 

(4.7b) 

2 1 % ^ 

6 -' ; a = Mg = T+a uJm ( To P of the acoustic branch), 

and a = Mg/J^ ■ (4.7d) 

The frequencies o) a andh Q get interchanged when M 1 > Mg . 
llo’-o the two branches are each three fold degenerate. The 
Green’s functions for such crystals become wholly real inside 
the gap besides the out of band regions and complex within 
both the allowed bands. 



The perturbation and Green’s function matrix for the 
diatomic simple cubic lattice can be treated now on the sane 
footing as that of the simple cubic lattice after the lit' 
transformation. Since any physical property like that of 
infrared absorption depends on the unmodified Green’s funct- 
ions, one has to transform them properly before making any 
quantitative estimation. With the help of the transformation 
equations(4.2o) and (4. 2d) one can construct matrices which 
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transform the matrix Mu: 2 to Mcj® as 

(M.^. 2 - A) = S(£^0 2 -• 1) S" 1 (4.8) 

* */ 

wlioi o S and S are diagonal matrices whose elements are 
given by 


and 


6 i,J /C( “l‘ t,2 - 6Y)l5 3,e T en + ( M 2 uS - 6 V ^.oda^ ’ 

(4.9a) 

S i,jT = )*« 3 , e W<%« 2 -SV )*6 , , oaa ] . 


j,odd J 
(4.9b) 

Honce the Green’s function matrix occuring in equation 
(2.49) of power absorption is expressed as 


(Mc<J 2 -A) 1 = S 1 (LIu> 2 -A) 1 S , 


or 


A ^ *1 yC A ' 

G = S 1 G S 


(4.10a) 


(4.10b) 


In case of charged impurity only one element of the pert- 
urbed Green’s function is needed for calculating the 
absorption co-efficient. Here the entire lattice consists 
of ions and the sum over j and j/ should be carried out over 
nil the lattice sites. The charge of any ion is 


e. - (>l).?e * 
A 


(4.11) 


So the absorption line shape expression becomes 


BM=fc e s E a «J C - 


6, 


3 1 0dd 


•* 


.(M^ 2 -6# (M^> 2 -6y) 
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[(M 2 ,JM Y )^ evm -<l^ 2 -6Y> 


V,°m ] • < 4 - 12 > 


The method of partitioning as discussed in Chapter III 
helps us to sun over j and j/ocnveniently. Here the 
equation (3.13) is rewritten in a slightly different form 
in order that the unperturbed Green’s function separates 
out. 


G 



A & A ✓ 

+ H S 


(4.13a) 


with 

/-g V ( I+g* fr 1 o 0 .. \ 

He/ 

-4i p ' ( I+g * p) -1 ° 0 • • i 

V . * • • • j 

/ 


(4.13b) 


Now substituting the first term of equation(4.13a) in the 
expression for power absorption and making use of the 
explicit form of the Green’s function, the summation over 
all the lattice sites can be done easily. The contribution 
of perfect crystal to infrared absorption has the following 


form 


B p M 


- Mn +M r> 

\ ^ 2 ( i or ) 1111 


Lt}‘ 


ie 


% 

X 


M-, +M 9 


6 < ) 


(4.14) 


The 6-function shape of the absorption is the well 
known phenomenon in ionic crystals. In case of perfect 
crystals there appse^s a continum of weak absorption with 
a strong well defined maximum at the transverse optic mode 
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frequency. The continum arises due to anharmonicity terms 
or the higher order dipole moment of the crystal. But 
these processes involves more than one phonon and bring 
temperature dependence to the absorption. Both these 
processes are similar except one fundamental difference. 

Although a continuous absorption spectrum is obtained due 
to higher order dipole moment the infinitely strong absor- 
ption or the 6-function behaviour is reduced only by the 
anharmonicity which couples the phonon giving a finite life 
time of this mode to decay into the neighbouring modes. 

On the otherhand the defect also relaxes the 
selection rules producing absorption at all frequencies. 

In the harmonic approximation the dispersion oscillator 
predominates in its neighbourhood. Any defect- activated 
mode arising in that region will be dixxiowat to distin- 
guish. Only th^so modes which lie far away from the 
Reststrahl region can be studied properly. This makes the 
impurity resonant modes in the acoustic branch much more 
amenable to experimental study than those inside the optio 

branch. 

Ihe second term of equation( 4.13a) is the contribu- 
tion of the impurity to infrared absorption. mien^HG is 
substituted in equation(4.12) the summation over /is carr- 
ied out easily, beoause it occurs in the exponential m the 
integrand of the integrals for the element of G°. Ihe 
structure of H matrix as given by equation(4.13b) is suoh 
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that it has non vanishing elements only upto column whose 
index corresponds to the last lattice site affected hy the 
impurity. This enables uo to write the impurity absorption 


as 


B I (eJ) - i e 2 l 2 t o Im jg [ 


! )qj 2 -12y+2a 


a - 6y 




, , (%+M 2 )4j 2 -12Y-2a ^ 


2E 


+ 


a + 6y 


l— t H-? n' 




z 


a + 6y • 2,j±' 


•“■‘i,- * a.*s ,] 

a— 6y ■“*»- 2.* 


(4.15) 


where p runs over all the affected lattice sites and a is 
given by 


a = [(M 1 ^ 2 -6y)(H 2 to2-.6y)]^ 


(4.16) 


The sum over j con he performed hy malting use of 

i 

the partitioned form of 1 1. For example the first term in 
the above equation becomes 

x e i(n+p+p)7t .(4.17) 

where j. = n + £ and n runs over all the centers of the 

partitioned blocks of the Green's function matrix. By using 

x- 

thc explicit formula for the element of the sum is carried 
over n+]i and the result is obtained in tne following form 


1>& 


J i(.‘i+^)% f 
i 



(4hn' [ ^ 
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hzZL.lv (i+CT)" 1 ],, 


a-by 




£»Ji 
(4.18) 


Tlie other terns of equation(4.15) can he simplified in the 
similar manner to give 


_ , v 122 1 r 12y+2a sr~ > 

B t (co) = e 2 E 2 cJ Im -r± [- A ... 2 ! > 

4a ( a-6y ) 2 




lviUsh-% , e itt(w) + (Hx^z^-izr-aa ^ 

( a+6y) 2 


-j— M-i )t^ 


[id+g P) ^'4- - f , 4^, C® < I+ s s> V 


a -36y 2 


M 


(e i7t ^ + e i7 ^)J . (4.19) 


Nov; the sum extends only to those lattice sites which are 

affected hy the impurity and the total number ‘is -limited 

in this case to seven. Instead of summing directly over 

ji and j /we adopt the following procedure so that the result 

can he expressed in a more convenient form. The symmetry 

of the defect lattice and the V matrix which block diagonal- 

isos the submatrix (I+g*p) were discussed in Chapter III. 

Here the same V matrix can block diagonalize the matrix 

s 

p* ( i+g'p r 1 into its various irreducible parts. Then the 
summation over £ and £ is carried out for each partial wave 
blocks. The various sums that oocur in the above equation 


N, 


are 



2Z 

*2L, Cp* (i+s* pV 1 ] 


• 



• 

• 

CD 

O 


i 

D S 

(R-Q+6S), 

(4.20a) 

i±>i = 

1 

5 s 

(R+Q+6S), 

(4.20b) 


2 

CpWp*)- 1 ] '(s-ss). 

■ii.Ji £jJ! ' D 5 

( 4 . 20o ) 

whore 


• 

R - (^[a+ AM-jU) 2 (Mjcu z -6y)g^] , 

(4.21a) 


Q = 12 Ay (1+ AMjCi ) 2 ) , 

{4.21b) 


S = - ^ ( 1+ & MjU) 2 ) . 

(4.21c) 

and 




dS=1+ ^ ^Cl^) -V ♦ Btf+gj. 4?)- 


M-,o3 2 

C 1+ T;>W 1- 

(4. 21d) 

Here 

*g i 

D is the s-wave part of the determinant 

I+g p . The 


sum over the other two types of modes vanishes indentically. 
So the infrared absorption takes place only through the 
totally symmetric mode in ionic crystals and is similar to 
the other cases discussed in Chapter III. The physical 
picture for such phenomena is as follows* * 


The p and d-like modes involve the motions of the 
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nearest neighbours of the impurity. For ionic crystals 
they oarry similar charges hut move in such a way that the 
resulting dipole moment as go oi at ed with these modes vani- 
shes. For example in p-ljJke motion (Fig. 2b) a pair of 
ions on both sides of the impurity move in opposite direct- 
ions by the same amplitude cancelling each other's dipole 
moment. 

In d-type motion (Fig. 2c) either all the six or 
four of the nearest neighbours take part in the vibration. 
For the first case four of the atoms move in one direction 
and bho remaining two move in opposite directions with 

* 

amplitudes twice as that of the others. In second case 
two of the atoms move in one way forming a pair and the 
other pair moves in the opposite direction. Ihe resulting 
dipole moment in each of those cases become zero. Hence 
these modes cannot bo excited by the infrared even in the 
ionic crystals. Fhe possibility of exciting other modes 
by the infrared will ariso if the charge symmetry in the 
surrounding atoms near the impurity is destroyed. The 
process of direct calculation cannot reveal the physical 
situation as clearly as discussed above. 

So finally the absorption formula reduces to 

BjM = £ e 2 E 2 <0 Im . (4.22) 


where 
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N 



£M x 

-a+crH -Jp 




u ' 2 

— — — y 


-«4)+a(l+<r) (-^) (^) )l^3 . 


m 


(2.23) 

In deducing the above expression the following Green's 
function relationships (similar to those of the simple ouhio 
monatomic lattice) have been used. 

ag 0 + 6yg x = 1 , (2.24a) 

s*l + t(Sq + 4 §il + §2 ) = °- (2.24b). 


' ' 4.3 Density Of States 

The squared distribution function of frequency is 

given by 

/-(v 2 ) = W Tr (MG) = m Tr (^G) ; (2.25) 

where l' = S M S’" 1 . (2.26) 

Equation(2.25) can be written as 

- IrM°&°P (l+s%) -1 e°] .(2.27) 

whore ^ ( v0 2 ) is the density of states for the unperturbed 
lattice -and is given by 
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P(( M S ) = «2L fp r ffr°pOx _ im , M y 

It : Nic • Lr( ' L G ) - 2^ (l^P + -p^)g 0 (2.28) 

Hence the change in the density of states is given by 


' ) = m C hMipG(0,0) - M°G°P (I+gV)" 1 ^] 


* x 


_ In 


M, 


J T I r ~ w ' K a- _ 

Brf G(0,0)+J(pM 1 + ~|)PrG°p i< (I+G 0 P)“ 1 

“i 

- ^)Tr R P ( I+G 0 ?)- 1 ]. (2.29) 

whoro wo hnve used the following relationships. 

„ r M 0 M_ 

M i,i , = 4 a i,i /C(fiM x + r> + (PM 1 - r> • **] • 

(2.30a) 

2 -/ G °CL»£> G °(£>i) = - h G °Q>i) = - G °(i,i) , (2.30b) 


and 


RQ*JL) = 2j G °(i»sJ G °(^i) e lnil = ^ G°(i,2)[(-l)3 


11 


+ (-l) 3 ] . (2.30c) 


The traces of the matrices occuring in equation(2,29) can 
bo evaluated in al a following way. The second trace is 
easily calculated by the use of the partitioned form of 
their respective matrices which reduce the dimensionality 
of the resultant matrix i.e. 


Ir if j? V ( I + 0° P)" 1 = / p* (I + g v )~ X - (2.31) 


* ** \-*l 
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wht.i c r is the first hlook of the partitioned R* matrix. 
Then one oan calculate the trace with the help of the group 
theoretic block diagonalization procedure. The other trace 
is expressed as 


Tr G°P (I + G°P r 1 = Tr log (I + G°P* ) - Tr 


O A 


*^r\ a 




da 


= T r f a lo g (I + &Y) - [ 


(M^ 2 ~6y) 


/iY 6 y(M 2 -M.) * M 9 

+ Tri . J x G ' r> 


a (M^> 2 -6 y)' 

^ + ® ( Ji » jti )»(2.32) 

-fov) r II 


|1[ 


M, 


(M 2 u^-6y) 


M- 


w’n to 


& (j±*j±) = Tr (I + g*p*) g* - , 


(2.33a) 


and 


0 (0,0) = (g* Q + ^- 4l g£)/ D s .. (2.33b) 


How substituting equation(2.31) and equation(2.32) 
in equation^ 2. 29) and by making use of the explicit forms 
for equations( 2.33, 2.31), the ohange in the distribution 
f ■unction for the squared frequencies is expressed in a 
convenient form as 


= H Is log Bet (I + ® ® * 


/ 

40* 
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^ Im r D 


j.ia r JJ „ tf 
Hit *• - + 3 — + 2 


(3.34) 


whore 


I + g p 


* f 3 ;2 

D S pP D d , 


(2.55a) 


^ ■ 1 - («; - 4 ) , 

3)4 = 1 - ^ (e 0 + 4 - 2 «h ) 


(2.35b) 


(2.35c) 


i. 4.4 -Discussion 

I'be Green’s functions for the diatomic simple cubic 

lattice can be expressed in terms of the Green’s functions 

of the monatomic simple cubic lattice for which extensive 

tables are available and hence, the absorption line shape 

function o .an be plotted for various values of impurity 

*- 

parameters. The denominator h s is real inside the forbidden 

gap besides the out of band region and the poles lying in 

those regions give the gap and local modes respectively. 

* 

This local mode frequency when arises due to hydrogen impu- 
rity in alkali halide crystals is known as IT-center and 
wore observed in the optical absorption experiments as 
they are infrared active. ■ 

Inside the band D s is complex and the poles of its 
real part may lie either inside the acoustic or inside the 
optio branch depending -on the perturbation parameters 
( 9 21) and the masses of the atoms constituting the 
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1 at tide. As regards the aooustio branch is concerned one 
gets the usual, low frequency resonances for heavy mass def- 
ects or for impurity atoms associated with softeni ng of 
nearest neighbour force constants. This behaviour is exactly 
similar to that of monatomic lattices. But inside the optic 
branch things happen in a altogether different way depending 
on the masses and M2. For and the substitutional 

impurity replacing one gets resonances when the mass of 
the impurity is less than that of the host atom. As the 
masn of tho impurity decreases this resonance frequency 
shifts towards the top of the optic band and a further 
decrease in the mass shifts the mode out of the band and the 
localized mode is formed. Here the impurity host nearest - 
neighbour force constant plays the usual role i.e. decrease 
of force constant lowers the resonance frequency. For M-^Mg 
tho light dofoct at site produces resonance which moves 
towards tho bottom of the optic branch with decrease of 
mass of bhe impurity. The effect of force constant is also 
reverse here because an increase of nearest neighbour force 
constant lowers the resonance frequency. In case of reso- 
nances in optic branch, one faces the difficulty to study 
the line shapes theoretically as the Eeststrahln absorption 
predominates and this has been pointed out earlier.- However, 
one oah estimate the nearest neighbour interaction from the 
infrared measurements which are either in the acoustic or 
in the outofband region of the frequency spectrum. 



..87 


Ihc first experimental evidence of impurity acti- 
vated mba-'d resonant modes in mixed alkali halides were 
obnotvod by Slovers 6j . He has observed a strong absorption 
lino in bhe fax infrared region (at 33.5cm” 1 ) of the 
electromagnetic spectrum in K2r:AgBr crystals (silver 
concentration is of the order of IQ 18 per cm 3 ). Similar 
observations were made by him for KCl:AgCl and KIsAgl. His 
subsequent study 69 of different impurities in KBr was 
striking in the sense that he could detect a very low 
frequency impurity mode due to light impurities like Li 6 * 7 . 

Simple mans defect approximation in all these cases 
cannot explain the results. Sievers and Takeno analysed the 
KBrsLiBr spectra on the basis of a simple cubic monatomic 
lattice with the defect as the only charged atom in the 
crystal and employing Debye approximation for the evaluation 
of the Green’s function. They could obtained a fit by vary- 
ing the force constant parameter. It required an enormous 
reduction of impurity host force constant to obtain the 
observed resonance frequency in the low frequency region. 

Haro wo have fitted our results to their experimental data 

6 7 

for KBr orystal with Li * ' and Ag impurities. Ihe absorpt-. 

6 7 

ion line shapes for KBrsLi ’ ' and KBr:Ag are plotted in 
Pig. 7 . 30 able Ho, 1 gives the experimental data on the 

resonance frequencies and the values obtained for the force 
constant parameter (^). 30he results on the basis of a 
simple cubic lattice with a charged impurity are also given 
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Resuls for various infrared impurity lines in KBr 
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Other relevant data for potassium "bromide crystal. 

(i) Mass of potassium ion = 59*1 amu 

(ii) Mass of bromine ion = 79.9 amu 

(iii) cr = M 2 /M 1 ~ 2 

(iv' Wavelength of maximum absorption = 88.3 p. 

The results for Li^*^ are in agreement with that of 
Sievers and Takeno who had utilized the Debye model. Since 
the electromagnetic radiation excites the impurity? Inb and 
modes in these systems at long wavelengths the details of 
the structure of the lattice become less important. For 
such frequencies the difference of masses of both the ions 
seems insignificant. This is primarily the reason one gets 
almost identical results for the force constant parameters 
in KBr:Li by treating the lattice as monatomic simple cubic 
But there is a noticeable difference of the two models for 
KBrsAg. This makes it possible to use a Debye model calcu- 
lation for low f??equency resonances. For higher frequency 
this model becomes inaccurate. 

From the above analysis of experimental and theore- 
tical situations of infrared absorption, in mixed alkali 
halides certain inferences concerning the impurity host 
interaction can be made. There is a remarkable softening 

of impurity host interaction for small ions "like Bi + in~KBr 

T 21 2 2 

and KOI. Thermal conductivity- 1 -^- measurements at low 
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temperatures show dips indicating resonance scattering of 
phonons at very low frequencies. The specific case of Iii + 
impurities in KOI and KBr were widely studied in order to 
determine the impurity configuration in the host lattice. 

Besides the thermal conductivity and infrared 
absorption other experimental observations like that of the 
electro caloric and dielectric measurements indicate a 
peculiar type of lattice distortion that arises when the 
size of the impurity ion is much smaller compared to the 
size of bbo host ion it replaces. Particularly this is so 
in case of Li* impurities in potassium halides where the 
volume of the potassium cavity is nearly 8 times larger than 
that of the lithium (the radius of ! + = 1.3 3 A 0 , the radius 
of Li = .68A°). These experiments suggest that instead of 
the normal lattice site as the equilibrium position, the Li + 
has a number of equivalent off site p^sibxons and it tunnels 
between these equilibrium positions. In this respect the 
behaviour of Li + ion is similar to that of the molecular 
impurities (GW*", H0 2 “, 0H~) when substituted in the crystal. 
The dipole moment, polarizability in case of monatomic 
impurities arise due to its off center positions whereas 
the molecular impurities possess them intrinsically due to 
their internal modes. 

The earliest theoretical investigation concerning 

125 

this central instability was carried out by Mathew * He 
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lias calculated the energy of the impurity hy taking into 
account the coulomb-, repulsive and electronic polarization 
terms and observed minima for Li + in KOI along <<100/ direc- 
tions. Dienes et al"* -28 performed a more elaborate calcu- 
lation which supported the multiwell shape for the potential 
well. Bowen, Gomez- and Krumhonsl 127 * 128 have extended it 
further by assuming a three dimensional model potential 
having the octahedral symmetry which may have 6 wells along 
<100} , 8 well along \111> and 12 wells along <110> direct- 
ions. By making use of the ICAO method and group theory 
they have calculated the ground state splitting due to the 
tunnelling motions and some of the experimental data were 
fitted to the various transitions. 

The multiwell shape of the potential inside the 
cavity is a valid description for small ions. For other 
impurities like that of Ag (ionic radius of Ag + = 1.26 A 0 ) 
the matching is more or less balanced and the usual descri- 
ption of the substitutional defect holds good. 



CHAPTER V 
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DEBYE- WALLER FACTOR 

In this Chapter we shall discuss Briefly the 
dynamics of an individual atom in connection with the 
Mossbauer effect. In most cases of resonant emession or 
absorption of y-rays Toy nuclei hound in a crystal the 
emitting or the absorbing nucleus happens to be an impur- 
ity embedded in the crystal and hence in analysing the 
results for such cases the effect of impurity on the lattice 
vibrations of the host crystal should be taken into consi- 
deration. 


The probability that a resonant y-ray emission 
or absorption Till take place without creation or destru- 
ction of a lattice vibrational quantum of energy is known 
as ’Debye-Waller factor* in analogy with the corresponding 
factor which determines the intensity of the scattering 
of x-rays or neutrons in crystals. This probability factor, 
its temperature variation and rhe temperature dependence 
of second order Doppler shift of the zero phonon line are 
directly related to the dynamics of the resonant nucleus. 
Since the Mossbauer active nucleus is a foreign atom in 
the crystal, its interaction with the host may be different 
from that between the host atoms and so the dynamics 
depends significantly in some cases on this parameter 
besides the mass difference. This makes it possible to 
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utilize Mossbauer effect as a tool for estimating impurity- 
host force constant. 


By general arguments of momentum conservation 
and translational invariance Lipkin^^ has shown that the 
relative probability of a Y" r ay emission or absorption 
which accompanied by a transition of the crystal state 
from j.E^.y to j is given by the matrix element 

P(E f ,E i ) =j(E f | ± yj S . (5*1) 


where '£ is the wave vector of the y-xay and R(jj.) is the 
instantaneous position vector of the resonant atom. Since 
it is assumed that the crystal state is unchanged during 
this process, the probability is determined by P(E jL ,E i ). 
But the initial state is quite arbitrary and so this 
probability should be averaged over an equilibrium distri- 
bution of initial states i.e. 


f = <P(E.,E i )> 



(5*2) 


The evaluation of this averaged quantity has been discussed 
by Van Hove. 1 ^ 0 For harmonic crystals this is directly 
related to the mean square displacement of the impurity as 


f * 





* 


(5*5) 


The exponent of the Debye-Waller factor gets 
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changed if the anharmonicity is taken into account. But the 
most dominant factor is still given hy the above expression. 
Evaluation of ^u 2 (jj.,0)^ has been achieved in various ways 
t airing into consideration the lattice models for host cry- 


stals. Host of the calculations were performed for Debye- 
solids with an effective Debye- temperature i.e. 


= 


_j 

4mkO, 


[i ♦ ( 1 2 


-fA 

- 5 -) ; 4 s - ] 

f e x -1 


(5.4) 


where the effective Debye- temperature is given by 


Q 


f “ 




(5-5) 


Here the impurity parameters enter into the calculation 
through the effective Debye-temperature* Erom the experi- 
mental value of Debye-Waller factor one can estimate the 

effective Debye temperature and hence the impurity host 

Idl 

nearest neighbour interaction. Bara et al have measured 
Debye-Waller factor for Ee^' in sixteen different host 
crystals and car.e to the conclusion tb „t the binding between 
the impurity and host remains the sane as that between the 
host atoms. Beoently hus.-sl nun"' , " et al have studied 
the s.ame problem in some of the above hosts like that of 
platinum, palladium, and copper for a wide range of tempera- 
tures and they could obtain a good fit by assuming the simp- 
lified model of an Einstein crystal with an effective Einst- 
ein temperature. The force constant was seen to increase by 
20 '/< in case of platinum and palladium and 20 '/- less in 
case of copper, 

This suggests the possibility of estimating the 
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impurity host force constant "by taking suitable models 
for the host crystals. She models in lattice dynamical 
calculations are difficult to handle as they demand compli- 
cated numerical work. Here we use the scalar models as 
has been discussed in earlier Chapters to fit up some of 
the available d?ta and estimate the change in the nearest 
neighbour inter action. Shis calculation, though based on 
an unphysical model, is presented here to throw some light 
on the nature of the computational problems arising in 
the process. 


She impurity is .assumed to be at the origin of 
the co-ordinate system. Using the time independent displ 
acement into correlation functions as given by equation 
(2.47b) the mean square displacement of the impurity is 


given by 

^U 2 (0,0)> = 


c <V) 

lim •»— tCr( 0,0, +iE)—G(C), ) 0,fi)“iE) J 

-jo 

<XJ 

~ f d'/> coth Im G(0,0,4)-ie) . (5.6) 


Shis is the general expression. She evaluation 
of the mean square displacement can be achieved in two 
ways. In first method we express the (0,0) element of the 
perturbed Green's function in terms of the unperturbed 
Green's function as has been discussed in Chapter HI. -fin 
additional term arises in this case if the denominator 
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vanishes outside the hand producing a local node. 
equation(5.6) reduces to 


x> 


/ a/« A x\ ooth Img 0 

/u 2 (0,0)>= — ■ ■ X . cU-» ^ 

X ' 71 J B^+Df^) 


+ 



Then the 


(5.7) 


where NG-Op) is the numerator of the (0,0) element of the 
perturbed Green’s function, cjn is the local mode frequency. 
The first term is the contribution of all the band modes 
and the second term is due to the local mode. At finite 
temperature the local mode contribution is not negligible. 

It has been shown^ ’ that in the presence of a local mode 
a correction term to the lebye-Waller factor is to be 
multiplied i..e. 


f = 


-K 


X' 


t 2 (o,o))> 


0 


(z) 


(5.8) 


where I 0 (Z) = Bessel function of imaginary argument and 

Z- 2®] (5.9) 

einh (S^/ 2 ) 

W(1M) is the contribution of the local mode to the exponent 
of the Debye-Waller factor at T = 0. At zero temperature 
the correction factor reduces to unity and at higher tempe- 
ratures it enhances the value of f. 
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2he Green's functions axe calculated by the 
method suggested by Mahanty.^-^ The Debye- char act eristic 
temperatures'^'* axe used for the host crystals to fit up 
the experimental data. 


The second method ho\vever is more suitable for 
such calculations as one need not require the computed 
values of the Green’s function for the entire band. Using 
the expression for the perturbed Green’s function (equation 
2 . 51 )i the mean square displacement is 



ooth(p.y2) j ( 2 )/(£>. 

( 5 . 10 ) 


(PfcV' 2) 


s 


etc; 

At high temperatures the expansion for ooth in a power 
series is well known and is given by 


oothCpuig/a) ^ ^ , 

(^^3/2) " (P - S ) 2 


• l l 
+ 3 - X5 


C^r 


< Tt] 


+ #. 

(5.11) 


Substituting equation(5.H) in equation( 5.10) and making 
use of the closure relation as given by equation( 2.40b) the 
mean square displacement comes out to be 


u 2 (0,0)> - 

\ 


-i— [ - m'g( 0,0*0) + T2 P 2 


d(.q..q) 

720 


P 4 



p 6 + .. ]. (5.12) 


Here only the (0,0) element of the perturbed Green’s function 
at one frequency (co =0) enters into the calculation. Other 
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quantities can oe calculated easily. In the scalar models 

the various quantities that occur in the above equation are 
given by 


1(0,0) = = 1L.1L ( 1 IL) 

“ “ W 2 1 y * 


(5.12a) 


D 2 (0,0) = i ( SI) 8 (n + 111) =1 (L^) 2 ( 1_ 2) 2 
’ ” n M' ^ M ' n 1 2 ' ^ M '* y' 


(n+ §-), 


M' y ; 

(5.12b) 


and 


'J(O.O) = g 0 + i ( I - I ( ) = _ I ( S + Ij . (5.120) 


The final expression for the mean square displacement becomes 


k„T 


/u 2 (0 o)^ = -21 Tx+ 1 + i- M-L^ 2 1_ 

< U» U L nY LX+ ,+ 24 mA m) - o* 


M__„ 

ny y/ ' 24 M' v T' “ 2880 M' 

(HL 1) /JM + i m_/M x\ z 

K M‘ • y' ' T' + 241920 y' # 

( n + j|) (~rf?^) 0 + ... ] . (5.13) 


TT 


where kg is the Boltzman constant. At high temperature this 
expression holds good and for many hosts the limit of tempe- 
rature ( T ^ ) is sufficiently low so that the Debye-Waller 
factor can be computed for a wide range of temperature. The 
value of x is determined from the (0,0) element of the 
unperturbed Green’s function at^= o. 


Platinum, palladium and copper hosts were studied 
by Nucsbaun et al. The Green’s function' for a fee structure 
in the scalar model are used to compute Debye-Waller factor 



.100 


and to estimate the ohonge in the nearest neighbour inter- 
action. For comparison a simple oubio lattice model is also 
taken as its values .are readily available. These two results 
together with the result of the direct expansion are tabula- 
ted in Tables 2,3,4 respectively for platinum, palladium and 
copper. 

In case of Pt and Pd the calculation agrees fairly 
well with the experiment even at high temperatures. The 
direct expansion method and the result by integration method 
agree reasonably at .almost all temperatures except at 
T = 20°IC (which lies in the region where the expansion is 
not valid) . In fee model the force constat change is 
around 70 /* whereas the simple cubic model it is about 50 /* 
However, there is a significant deviation for ou at high 
temperature. This discrepancy is still retained in the 
simple cubic structure. 

It is seen that in this model calculation there 
arises either a local mode or an inband resonant mode for 
the above mass difference and force constant change para- 
meters. Since the Debye-Waller factor is an integrated 
effect over all the modes bhe presence of resonances or 
localized modes in the spectrum would make it possible to 
use " the Einstein oscillator model as -has been taken by 
Hussbaum et al. 

In the absence of suitable lattice models no 

be made, but the results are of 


realistic estimation can 



(a) By method of integration 

(h) Dire ot' expansion method x = .2768 
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Table No 
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Table No 



i 



Table No 
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qualitative interest. The change of nearest neighbour 

interaction is significant when the mass difference between 

57 

the host and the impurity is large, like that of Fe in 
Pt and Pd. But for small mass difference the anharmonio 
interaction should be taken into account. 

The expression for the mean square velocity which 
determines the second order Doppler shift of the zero phonon 
line is obtained as 

<’4 (°.°)> “ 1 [ + T5 V 

* ^20 B aot (5 ‘ U) 

This expression is general and the estimation of 
nearest neighbour force constant can be made for any model. 
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CHAPTER VI 
CONCLUSION 

In this chapter the results obtained in the present 
investigation ruid their limitations are briefly discussed. 

The expressions derived for power 'absorption due 
to a substitutional impurity end the tonal integrated absor- 
ption ore general for harmonic crystal possessing a first 
order electric dipole moment. The absorption line-shapes 
arc evaluated theoretically in detail for charged impurities 
in otherwise uncharged host lattices (simple cubic, body- 
centerod cubic and face-centered cubic) and for mixed alkali 
halides (di.atomiG simple cubic lattice with alternate +ve and 
-ve charges) for scalar models, with the hope that the exact 
solutions in these ' idealised models would give some insight 
to the real problem. 

The group theoretic method clearly reveals the 
symmetry mode vibrations of the defect and its surrounding 
atoms which are affected. Out of all the symmetry modes of 
vibrations possible in the crystal only one mode namely the 
s-wnve vibration becomes infrared active. The inband reso- 
nances occuring in this branch have been fitted to the exper- 
iment ,ol data obtained from the infrared absorption in alkali 
halides. The tremendous softening of nearest neighbour force 

c n 

constant especially for small ions like Li * in KBr suggests 
tho departure of impurity potential from that of the harmonic 
approximation. As has be -mJLis cussed in chapter IV the 
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potential inside the cavity has a multiwell shape for small 
ions and quantum mechanical description of tunnelling motion 
•is reasonably adequate to describe its dynamics. However, 
the interaction of such an atom with the rest of the crystal 
should be taken into account for a better analysis. 

Using the experimental values for the Debye-Waller 

pr 7 

factor of Be 5 in Pt, Pd and Cu an estimate of the impurity- 
host force constant is made. 


All the results obtained are more or less qualita- 
tive, because any quantitative estimation should take into 
account the realistic features of the lattice i.e. the true 
density of stakes for the host crystal. The formalism used 
is very general. Further study can be made by using a 
polarisation model for the host lattice within the’ same 
theoretical frame work as presented here, but it involves 
far more elaborate calculations and numerical computations. 
Besides, the theory should incorporate anharmonicity, higher 
order electric moments and long range interactions between 
atoms for a complete description. The later interaction is 
of importance in case of ionic crystals. 


The electric field B of the electromagnetic radia- 
tion should be replaced by an effective field (or‘ local 
field) E3.0C inside ^ lie dielectric material, which takes into 
account the polarization of the medium. . This brings into 
the picture the dielectric constant of the crystal and the 
intensity of the source as parameters. 
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Appendix I 


The vibration of an impurity atom is related to 
the perturbed Green’s function of the lattice. As indicated 
in equation(3.16), the perturbed Green* s function 
0(0,0„...> 2 -ir.) which determines tbe vibrations of the impurity 
and hence the infrared absorption, equals [(I+g p)g]Q Q . 

The matrices g, p and hence (I+g p)g can be block-diagonal- 

A 

ized by a unitary matrix Y into submatrices belonging to 
different irreducible representations of the point group 
that occurs in the reducible representation generated by 
the affected lattice sites (equation 3.18), Denoting the 

A ^ A ~K\ K 

block -diagonal form of (I+g p)g by B, we have. 


B = Y + (I+g p) -1 g Y (A-l) 

B, 0 0 ,7\ fO 0 0./ 


M x i 

B 0 0 , . ! ! 0 0 0 .. 

£ j : 


• 

• 

O 

CM 

PQ 

O 

0 B, , # > # • . * , 

+ 

, . • • • 

* * . . , » . * •• 

' \ 


• • • • • 

1 

" h. J 


W / A ^ 




V 


From equation (A-l) we obtaih. 


(I+g p)g = Y B Y’ 


-2L 




A A A jl . 


(A-5) 


where each term on the right side of equation( A-3) belongs 
to a distinct irreducible representation. The (0,0) element 
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of G- now becomes 


G 


(0,0, <o 2 -i E ) = Z^V +) o,o • 

h 


(A-4) 


I'he V matrix is constructed as follows. The vector 
defined by, 


9 


/x(0) 
; KD 


(A-5) 


; X(n) 

\ 1 


I 


generates a reducible representation p whiGh can be reduced 
as. 


f= Z a r r ( ^ 


r 


(A-6) 


where a„ denotes the number of times the irreducible repre- 

‘ i — ’0 ^ *— • * 

sentation / occurs in f . If / is n - dimensional, the 
matrices for the operation 0 R (corresponding to the 

group element S) will have n columns and rows. Then the 

* A 

columns of the unitary matrix V, which are orthonormalized 
vectors, can be projected out of by the projection operatiif 
in the form 


i< M = 3 ± 


-Ji T>^(R)0 

g R ^ii wu R 


(A-7) 


where g Is the order of the group and the index runs from 1 

A 

to n . Thus the columns of V divide into sets of n.. vectors 

\l < r 

•* t * j \ . 
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for each and the members of the set are labelled by 

the index i. * 


In our notation the indices of the first element 
of the perturbation matrix P correspond to the impurity 
fc-ite. The elements of the first row of the V matrix which 
are given by ^'^(0) for various values of p and i can be 
evaluated trivially, noting that all operations 0 R of the 
group leave the first element X(0) of p unchanged. 
Therefore, 




(A-8) 


However, from the orthogonality relations for the matrix 
elements of different irreducible representations it follows 
that 

'4i )(H) = 0 * U-9) 

iv. 

for every p, except the totally symmetric one- dimensional 
representation, for which we get, 

I ^ (R)-, i = 1 . ( A-10) 

g R 

A 

Thus all the elements of the first row of Y are zero except 

A 

the first one, which is unity. Prom the unitarity of Y 
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matrix it follows that the elements of the first column 
(except the first element) are also zero. Hence we gel; 


Gr(0 


o,o,w*-fe) = 21j t5o,o 


= V(0,n) M u (n,m)Y + (m,0) 

(i,n, m ^ 


= Zj ^(0,0) ( A- 11 ) 

A 

The only M which has non-vanishing (0,0) - element is the 

A 

one that is generated by the columns in 7 matrix correspond- 
ing to the totally symmetric irreducible representation. 
Hence it follows that the block corresponding to the totally 
symmetric or s-wave representation contributes to the 
impurity vibration. 

A 

The explicit form of V matrices are given below 
(i) Simple cubic lattice 


A 



/ 1 

0 

0 

0 

0 

0 

° 

0 

a 

d 

0 

0 

0 

•® I 

; o 

a 

0 

d 

0 

c 

i 

“ e i 

i 

0 

a 

0 

0 

d 

-2c 

1 

/0| 

1 

0 

a 

-d 

0 

0 

G 

e l 

i 

! o 

a 

0 

-d 

0 

G 

-e . 

0 

a 

0 

0 

-d 

-2o 

oj 




2 ^ 


g 


x lu 


E. 


g 




(A-12) 



Ill 


(ii) Body-Centered cubic lattice 



(iii) Pace-centered cubic lattice 

1 0 0 0 0 0 0 

0 c Id 0 b Id h 

0 c b 0 -b b h. 

0 c b b 0 0 -2h 

• 0 c b -b 0 C -2b 

0 c 0 b -b -b h. 

V F =;0c0bb-bii 
lO c -b 0 -b b h 

| 0 c -b 0 b b h 

| 0 c -b -b 0 0 -2h 

0 c -b b 0 0 -2h 

1 0 c 0 -b b -b h 




( A-14 ) 
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Here 


a = ~— 

•u _ 1 

y D — 1 } 

0 = — 

v* 

r s 

1/12 


, e = — , 

h = i- 

V2 

2 

1T24 
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